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Abstract 

The  United  States  is  very  dependant  upon  the  use  of  space.  Any  threat  to  our 
ability  to  use  it  as  desired  deserves  significant  study  One  such  asymmetric  threat  is 
through  the  use  of  a  microsaiellite.  The  feasibility  of  using  a  microsaiellite  to  accomplish 
an  orbital  rendezvous  with  a  non-cooperative  target  is  being  evaluated.  This  study 
focused  on  identifying  and  further  exploring  the  technical  challenges  involved  in 
achieving  a  non-cooperative  rendezvous. 

A  systems  engi  neen  ng  analysis  and  review  of  past  research  quickly  led  to  a 
concentration  on  the  guidance,  navigation,  and  control  (GN&C)  elements  of  the 
microsateiliie  operation.  While  both  the  control  laws  and  orbit  determination  have  been 
previously  evaluated  as  feasible,  the  integration  of  the  two  remained  in  question.  This 
research  first  validated  past  efforts  prior  to  exploring  the  integration.  Impulsive  and 
continuous  thrust  control  methods,  and  linear  and  nonlinear  estimator  filters  were  all 
candidate  components  to  a  potential  system  solution. 

A  simple  yet  robust  solution  could  not  be  found  to  meet  reasonable  rendezvous 
criteria,  using  essentially  ofT-ihe-sheif  technology  and  algorithms.  Results  reveal  a 
simple  linear  filter  is  a  misapplication  and  will  not  at  all  work.  A  nonlinear  filter  coupled 
with  either  a  continuous  or  impulsive  thrust  controller  was  found  to  get  somewhat  close, 
but  never  close  enough  to  attach  to  the  target  satellite.  Successful  ON&C  subsystem 
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tntegraiion  could  only  be  achteveil  for  a  very  simple  case  ignonng  orbii  perturbauons 
such  as  tbe  earth's  oblaieness  A  top-level  system  architecture  for  a  non-cooperative 
rendezvous  nucrosaielliie  has  been  developeil  The  technical  complexity,  however, 
requires  more  complex  algorithms  to  solve  the  rendezvous  problem. 
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I.  loiroduciloo 


Background 

We  are  enienng,  or  have  already  eniered,  ^  which  the  uae 

of  space  wii  ( even  ^uch  [vofoutvj  influence  on  humnn  affairs  ihai  no 
notion  will  be  fully  able  to  conirol  its  own  destiny  without  significant 
space  capabilities. 

-General  Robert  T.  Ilerres  (USAF) 

Vice  Chairman  of  the  JCS,  1988 

Although  ihe  quote  above  may  seem  a  bii  dated,  we  are  just  now  beginning  to 
fully  understand  ihe  ramifications  of  such  statements.  The  United  States  has  for  some 
time  led  the  world  in  the  use  of  space-based  resources  for  military  as  well  as  civilian 
operations  Our  large  competitive  advantage  is  now  being  challenged,  however.  Gen 
Lance  Lord,  Commander  of  Air  Force  Space  Command  (AFSPC)  discussed  this  issue 
with  space  industry  leaders  in  November  2003.  The  Genera)  staled  that,  *'Our  adversanes 
-  and  even  future  adversaries  -  know  the  value  we  place  on  space  to  enhance,  improve 
and  transform  all  our  operations.  They  will  increasingly  try  to  deny  us  the  asymmetnc 
advantage  that  space  provides”  (Wi)son,  2004). 

As  the  potential  benefits  of  space  operations  are  more  broadly  understood,  more 
capabilities  are  being  transferred  to  the  ultimate  high  ground.  In  an  interview  last 
October  with  Inside  the  Pe}tia^on>  Lt  Gen  Dan  Leaf,  Vice  Commander  of  AFSPC, 
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describe>l}ijsi  rvow  de^^endeni  ihe  United  Slates  is  on  our  space  assets.  He  ctwse  lo  depict 
our  space  capabilities  as  '‘woven  inevincably  through  our  overall  nulitary  capabilities'' 
(Grossman,  2003al.  We  have  indeed  moved  much  beyond  ihe  first  space  war  of  Desert 
Storm  in  1991  when  OPS,  DSP,  and  exclusive  national  systems  were  used  to  support  only 
selected  air  and  ground  operations. 

Given  our  critical  dependence  upon  and  potential  threat  to  our  space  assets,  any 
focused  research  in  this  arena  may  prove  valuable  to  space  policy  maJters,  developers  and 
operators  alike  Gen  Lord  declared  that,  “(i  is  our  duty  to  preserve,  protect  and  defend 
the  high  ground  of  space  and  we  must  have  the  ways  and  means  of  detecting, 
characterizing,  reporting  and  responding  to  attacks  in  the  medium  of  space''  (Wilson, 
2(X>4)  This  research  efton  is  aimed  at  contributing  to  characterizing  a  speciRc  potential 
space  threat 

Space  Control 

The  concept  of  space  control  involves  both  offensive  and  defensive  activities  to 
ensure  a  desired  level  of  advantage.  From  the  very  indiscriminate  nuclear  systems  to  the 
Laser-focused  Star  Wars  initiative,  history  provides  a  colorful  review  of  space  control 
anempis.  Only  one  select  example  will  be  quickly  discussed  here.  Current  space  conuoil 
doctrine  will  be  presented  next  Finally,  the  above  will  be  used  to  put  this  research  into 
the  larger  context  of  current  Air  Force  counterspace  activities  This,  in  itself,  is  a  systems 
engineering  activity  as  a  key  step  in  the  evaluation  of  any  potenual  system  is  lo  view  both 
internal  and  external  environments. 
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Tradiiionally,  an  Anii-Saielliie  (ASAT)  syitem  hjs  been  viewed  a&  having  ihe 
purpose  of  negating  ihe  funciionjl  nussion  of  the  target  space  asset  This  can  be 
accomphsheij  by  various  tnettv>d«  Directing  eneegy  on  the  satellite  ftom  a  ground  or 
space-based  illuminaiing  device;  placing  co-orbiting  "mines”  In  space  adjacent  to  the 
largei;  direct  ascent;  achieving  a  co-orbii  with  the  target  satellite  and  "catching  up"  to  it; 
and  launching  a  device  from  a  high-akirude  aircraft  have  all  been  attempted  (Johnson- 
Freese,  2000). 

The  United  States  has  undergone  the  most  extensive  ASAT  development  activity 
Project  SAfNT  (SAtelhte  INTerceptor)  began  in  the  late  l9S0's.  The  program 
extensively  covered  a  wide  range  of  technologies  for  interception,  inspection,  and 
destruction  of  enemy  spacecraft  (SAfNT,  2003 ).  The  Concept  of  Operations,  or 
CONORS,  entailed  rendezvous  with  a  target  satellite,  inspection  with  television  cameras, 
and  then  disabling  ii  somehow.  Project  SAINT  was  restructured  several  times  and 
eventually  canceled  in  1962  before  reaching  operational  status. 

The  doctrine  of  Space  Control  has  only  recently  emerged  within  both  the  DoD 
and  the  Air  Force.  DoD  3100.10  defines  Space  Control  as;  "Combat  and  combat  support 
operations  to  ensure  freedom  of  action  in  space  for  the  United  States  and  its  allies  and, 
when  directed,  deny  an  adversary  freedom  of  action  in  space  "  It  further  delineates  Space 
Control  mission  areas  to  include  surveillance  of  space;  protection  of  U.S.  and  fnendly 
space  systems;  prevention  of  an  adversary's  ability  to  use  space  systems  and  services  for 
purposes  hostile  to  U.S.  national  security  interests;  negation  of  space  systems  and 
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services  u&ed  for  purfo&es  hostile  lo  U  S.  national  security  interests;  and  directly 
supporting  bsiile  management,  command,  control,  communications,  and  intelligence. 

The  Air  Force  uses  Air  Force  Doctrine  Document  (AFDD)  2-2,  "Space 
Operations,”  to  specily  the  approved  methods  and  means  of  conducting  counterspace 
activities.  The  function  of  Counterspace  is  assigned  to  fulflil  the  Space  Control  mission 
area.  AFDD  2-2  describes  Counterspace  Operations  consisting  of  those  operations 
conducted  to  attain  and  maintain  a  desired  degree  of  space  superiority  by  allowing 
friendly  forces  to  exploit  space  capabilities  while  negating  an  adversary's  ability  to  do  the 
same. 

As  with  most  functional  areas,  Counterspace  is  divided  into  otFensive  and 
defensive  components  Offensive  Counterspace  (DCS)  operations  preclude  an  adversary 
from  exploiting  space  to  his  advantage  (AFDD  2-2,  2001 ).  The  usual  contmuum  of 
Deception,  Disruption,  Denial,  Degradation,  and  Destruction  are  available  means. 

Defensive  Counterspace  (DCS)  operations  preserve  U.S  /allied  ability  to  exploit 
space  to  Its  advantage  via  active  and  passive  actions  to  protect  fnendly  space-related 
capabilities  from  enemy  attack  or  interference  (AFDD  2-2,  2001).  Both  active  (e.g. 
detect,  track,  and  identify)  and  passive  (e.g.  survivability)  techniques  are  promoted. 

Space  Control  cannot  be  effectively  achieved  without  both  robust  DCS  and  DCS 
capabilities. 

This  research  directly  supports  the  ability  to  conduct  DCS  operations.  Analyzing 
the  technical  challenges  arising  from  a  microsatellite  rendezvous  concept  helps 
characterize  feasible  adversary  CCS  capabiltiies,  and  thus  necessary  U.S.  defenses.  Lt 
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Gen  Leaf  recenily  discuased  a  future  a^iace  activity  identification  capability  gap  resulting 

a  multi-service  review.  Tlie  general  asserted, 

We  must  have  good,  dmely  ^ce  situational  (SA)  awareness  •  not  just 
because  of  our  i  ncreased  reliance  on  space  capabi  lity  and  the  complexity 
of  all  that  occurs  in  space,  but  also  because  of  potential  threats  to  those 
capabilities.  The  number  of  nations  that  utilize  space-based  capabilities 
and  the  way  that  they  are  used  are  both  expanding  So  we  have  to  ensure 
that  our  space  S A... doesn't  simply  track  objects  but  is  able,  in  a  timely 
manner,  to  recognize  changing  situations  in  space,  just  as  we  do  in  the 
atmosphere  or  in  the  sea.  (Grossman,  2003b) 

This  theoreiicai  capability  gap,  between  what  we  need  and  what  we  have,  is  not 
well  understood.  The  fact  that  we  do  not  have  a  good  characterization  of  feasible 
adversary  DCS  capabilities  leads  to  a  poor  understanding  of  the  gap  Better 
comprehension  ot'chis  potential  threat  is  the  aim  of  this  study. 


Problem  .StoteniendResctrch  Objectives 

The  overall  objective  of  this  research  efton  is  the  analysis  of  potential 
counterspace  threats  from  foreign  countries  or  organizations.  Coimterspace  operations 
that  are  possible  with  readily  available  technology  and  information  will  be  evaluated. 
This  efton  is  a  systems  design  study  on  a  potential  foreign  offensive  counterspace 
satellite  lo  identify  the  technical  challenges  arising  dom  rendezvous  with  a  non- 
cooperative  satellite 

The  specific  objective  is  to  determine  if  it  is  possible  to  design,  build,  and  operate 
an  ottensive  microsaielllte  using  off-the-shelf  technology  and  information  that  Is  publicly 
available.  The  microsaielllte  must  be  able  to  maneuver  to  rendezvous  with  a  target 
satellite,  mamiain  proximity  with  the  target,  and  perform  its  mission. 
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For  this  propel,  it  is  assumed  ihe  micrc«aie!li(e  wiM  be  placed  into  an  orbti  similar 
lo  dial  of  die  target  satellite,  approximately  1 000  km  behind  ii  in  the  aame  orbital  plane. 
The  microsaiellite  then  performs  rendezvous  maneuvers  to  approach  the  target. 

It  IS  further  assumed  die  microsatelliie  has  perfect  knowledge  of  its  own  position 
and  velocity  but  must  estimate  that  of  the  target.  The  microsatelhte  would  likely  begin 
with  an  orbit  solution  derived  ftom  off-board  sensors  As  die  microsatellite  approaches 
the  target,  on-board  sensors  would  detect  die  target  satellite  and  an  updated  orbit  solution 
would  be  calculated.  This  would  allow  the  microsatelhte  to  complete  the  rendezvous 
without  any  feedback  from  the  target  satellite. 

The  unique  aspect  of  this  problem  involves  die  use  of  an  integrated  estimator  and 
conuoller  to  more  closely  model  reality.  This  research  then  takes  an  additional  step  to 
make  a  systems-level  feasibility  assessment  of  the  proposed  microsatelhte  threat. 

Methodology 

A  high-level  systems  view  was  coupled  with  a  more  detailed  technical  assessment 
to  form  the  approach  to  answer  the  research  objective.  Systems  engineenng  design  tools 
were  first  used  to  identify  die  driving  technical  areas  to  focus  on.  Once  identified,  these 
guidance,  navigation  &  control  (GN&C)  algorithms  were  studied  extensively  to 
appreciate  die  evident  as  well  as  subtle  application  challenges  involved.  Finally,  the 
systems  view  was  again  taken  to  make  die  concluding  feasibility  assessment. 
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1(.  Literal  ure  Revle*' 


Cbtpter  Overview 

TtK  purfo&e  of  (his  chapter  is  lo  provide  a  review  of  microsaielliie  lechttology 
and  recent  rendezvous-related  research  efforts  Applicable  industry,  academic  and 
military  microsaiellite  eftons  are  presented.  Significant  challenges  regarding  ON&C  are 
highlighted. 

Previous  microsatellite  rendezvous  research  results  are  reviewed  nevL  The  focus 
is  on  AFfT  work  leading  up  to  this  research,  and  supplemented  were  appropriate  The 
conucl  taws  and  orbit  determinaiioivnavigaiion  rei^utred  to  support  a  non-cooperative 
rendezvous  make  up  the  primary  body  of  research  drawn  upon  in  developing  (he  starting 
point  for  (his  pro^ct. 

Literature  reviewed  was  primarily  limited  to  open-source  as  (he  feasibility  ofa 
relaiivety  low-tech  solution  using  off-the-shelf  technology  and  publicly  available 
mtbrmaiion  is  being  evaluated.  The  review  highlights  several  key  findings  as 
summarized  below  The  microsatellite  industry  is  rapidly  becoming  capable  of  providing 
system  solutions  to  a  very  diverse  set  of  problems,  to  include  space  control  applications. 
Rendezvous  with  a  non-cooperative  target  is  a  non-trivial  operation,  but  key  elements  of 
conutil  and  orbit  determination  have  been  separately  demonstrated  Little  research  has 
been  accomplished  to  evaluate  systems-level  feasibility  from  a  systems  engineering 
perspective. 
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Vflcrosalelliies 


A  imcrosaielliie  is  a  small  soidlne  generally  considered  lo  have  mass  less  ilian 
100  kg  They  are  rypically  more  economical  lo  develop  and  operate,  and  quicker  from 
concepi  10  operation  compared  lo  tradiuona)  satellites.  There  has  been  considerable 
etTon  m  ihe  research  and  development  of  capabilities  using  microsaielliies  of  laie. 
Common  uses  include  visible  sensing,  mulii-specu^  imaging,  radar,  infrared, 
commimicaiions,  and  navigation. 

Surrey  Satellite  Technology  Ltd.  (S STL)  is  a  world  leader  in  the  develOfHitent  of 
microsaieiliie  technologies  Since  SSTL  spun  otT  from  the  University  of  Surrey 
Engineering  Department  in  191)5,  they  have  launched  roughly  one  spacecraft  per  year, 
pushing  small  satellite  technology  (Morring,  2003).  SSTL  claims  they  were  the  first 
professional  organization  to  offer  low-cost  small  satellites  with  rapid  response  employing 
advanced  terrestrial  technologies  They  indeed  have  an  impressive  track  record  in  an 
emergent  field. 

SSTL's  A(SAT-I .  developed  tor  the  international  Disaster  Monitoring 
Constellation,  has  successfully  completed  over  one  year  of  operations.  Imagery  derived 
has  been  useful  to  auihortiies  with  areas  of  responsibility  from  hydrological  mapping  to 
the  threat  of  locust  plagues.  The  AISAT-I  microsatellite  is  pictured  m  Figure  1  below. 


Figure  1.  A(SAT-I  Microsuielliie  (SSTL.  2003) 

The  European  Space  Ageacy  (ESA)  ia  alio  a  aigniijcjni  player  uaing  atnall 
seulbiea  for  advanced  science  eiisaiona.  ESA's  Projeci  for  On-board  Autonomy  (^ba) 
»  using  a  mscrosaielbie  lo  flight  teat  on-orbit  operational  autonomy.  Proba-K  iaunched 
in  Ociober  2001,  relumed  high-resoluuon  images  of  Earth  and  conducied  various 
radiation  studies  (Momng,  2003).  Proba-2.  sciteduled  for  launch  in  2006,  wtU  study  the 
sun,  providing  early  warnings  of  solar  Dares.  Frederic  Teeton.  Proba  prefect  manager 
notes  that,  ''small  satellites  have  proven  their  worth  for  rapid  leeimg  of  spacecraft 
lechniitues  and  onboard  instniments  They  can  also  support  dedicated  missions  very 
erricienily"  (SpaccDai)y,  2003)  A  Proba  microsatelliic  is  depicted  below  in  Figure  2. 
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Figure  2.  Anisi  Lnpres&ioD  of  Probe  io  Orbit  (ESA,  2003) 

The  Air  Force  Research  Laboratory  (AFRL)  ii  cuirenily  buiidiit^  and 
dcmunstrating  microsaieDiieiechnologioa.  SpeciflcaDy,  the  ExperimeniaJ  Spacecraf) 
System  (XSS)  Microsaielliic  Demonsiraiton  Project  me  (odes  two  very  applicable 
missions.  These  nussions  are  to  actively  evaluate  future  applications  of  microsatellite 
technologies  to  include:  uispeciion,  rendezvotji  and  docking;  repositioning:  and 
techniques  for  close-in  proximity  maneuvering  around  on-orbit  assets  (XSS- 10  Fact 
Sheet). 

XSS-IO,  Launched  in  January  of2(X)3,  commenced  on  autoitomoiLv  inspection 
sequence  around  the  second  rocket  stage,  transminmg  live  video  to  ground  stations.  Key 
technologies  demonstrated  mclude;  lightweight  propulsion:  guidance,  nat  igaiioa  and 
control  (GN&C);  and  integrated  camera  and  star  seiuor  (XSS-IO  Fact  Sheet).  XSS-10 
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achieved  ns  pnreoi^  mission  by  successfully  muicuvenny  frora  within  100  m  to  35  m  uf 
the  rocket  sta^e,  baLSung  away  and  repealing  (he  process  again. 

Work  00  the  fellow-oo  vehicle^  XSS-I I .  coniuues  XSS-1 1  is  to  further  advance 
technologies  and  lecluKfues  to  Increase  the  lexe!  of  onboard  auionomy.  One  of  the  major 
challenges  is  in  how  u>  sense  relative  position  and  velocity  when  in  proximny  to  another 
space  object  ( Paitch,  2003 ).  The  eftons  of  AFRL  con  firm  the  need  fur  complex 
navigeiiod  and  orbital  guidance  algorithms  onboard  the  spacecratl.  XSS-1 1  is  dJusu^icd 
in  Figure  3  courtesy  of  AFRL. 


Figure  3  XSS- 1 1  Operuting  a  Low-Power  Lidar  (Panch,  2003 ) 

XSS-I  I  IS  both  a  fast  paced.  JO-month,  and  highly  collaborauve  cfTon  The 
Space  Vehicles  Directorate  of  AFRL  is  paimenng  with  Lockheed  Martin  and  Jackson  & 
Tull  to  build  aod  iniegraie  the  microsaiellile.  (t  will  employ  a  sophisiicaied  three-axis 
siabilured  plaiform,  advanced  propulsion  system,  and  communicaiionx  subsystems 
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piishmg  the  scientirbc  envelope  This  will  all  lead  to  real-iime  streaming  video  of  the 
proximity  operations  being  sent  to  ground  operators. 

The  avionics  system  is  understandably  the  core  of  the  XSS-)  1  spacecrafL  The 
radiation  hard  Power  PC  750  processor,  develop  by  APRL  and  NASA,  enabling  the 
complex  data  processing,  guidance  algorithms,  and  onboard  autonomy  will  encounter  its 
first  tlighi  tesionXSS-l  l(Panch.  2003).  The  challenge  of  sensing  relative  position  and 
veloary  also  required  a  new  material  solution.  Due  to  the  lack  of  commimicaiion  with 
the  target  satellite,  AFRL  had  to  develop  alternative  approaches  for  relative  position 
determination.  The  active  sensing  system  selected  involves  a  high  tech  active  scanning 
lidar  ranging  system.  Complementing  the  active  system,  XSS-I I  will  employ  a 
combined  visible  camera  and  star  tracker  passive  remote  sensing  system.  Finally, 
onboard  iterative  trajectory  simulations  are  coupled  with  an  advanced  autonomous  event 
planner,  monitor,  and  forward-thinking  resource  manger  to  optimize  the  timing  of  rocket 
firings  (Partch,  2003).  The  development  and  miniatunzaoon  of  the  above  key 
components  required  significant  joint  research,  development,  and  integration. 

The  above  review  is  only  a  small  sample  of  current  and  projected  microsaiellite 
activity  in  industry,  academia,  and  military  arenas.  It  serves  to  support  the  argument  that 
microsatelliie  capabilities  will  continue  to  rapidly  increase  as  technical  hurdles  are 
overcome  Specific  to  the  non-cooperative  rendezvous  problem,  the  work  of  AFRL  is 
particularly  applicable.  ON&C  technology  maturation  is  an  area  of  intense  examination 
Advances  in  mmiaiurizaiion  and  the  proliferation  of  space  technologies  will  enable  many 
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less  knowledgeable  coimiriee  (o  contend,  as  unsolvable  probleins  of  today  will  be  taken 
for  granted  lomorrow. 

Rclevooi  Research 

A  wealth  of  previous  research  has  been  conducted  at  ART  regarding 
mtcrosaielliie  rendezvous  and  dockmg  operations.  Just  within  the  last  two  years,  work 
related  to  the  selection  of  tracking  and  orbit  determination  architectures;  rendezvous 
conutil  algorithm  development;  and  target  satellite  dynamics  modeling  for  microsatelliie 
docking  detection  has  been  accomplished.  A  recent  American  Aatronautical  Society 
paper  outlining  a  conceptual  design  for  the  GN&C  system  for  a  maintenance  and  repair 
spacecraft  complements  the  above  work. 

Cootrol  Laws 

Troy  Tschirhan,  a  former  ART  master's  student,  studied  the  control  bws 
necessary  for  achieving  rendezvous  with  a  non-cooperative  target  while  minimizing  fuel 
requirements.  The  relative  motion  of  a  microsatelliie  and  target  satellite  were  described 
using  Hill’s  equations  and  two  different  conuxiller  methodologies  were  investigated.  An 
impulsive  thrust  controller  based  on  the  Clohessey- Wiltshire  solution  was  found  to  use 
little  fuel,  but  was  not  very  robust.  A  continuous  thrust  conuxiller  using  a  Linear 
Quadratic  Regulator  (LQR)  was  found  lo  be  more  robust,  but  used  much  more  fuel.  The 
algorithm  developed  for  this  control  method  is  depicted  m  Figure  4  below. 
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Figure  4.  Linear  Quadratic  Regulator  Propagaiion  Algonthm  (Tschirharr,  2003) 


As  e  final  soluoon,  a  hybrid  controller  was  evaluated  which  uses  ihe  low  thrust 
Clohessey-Wiltshire  a^iproach  lo  cover  most  of  ihe  necessary  distance,  and  then  switches 
lo  the  Linear  <^adraijc  Regulator  method  for  the  final  rendezvous  solution.  Results 
show  that  this  apfiroach  achieves  rendezvous  with  a  reasonable  amount  of  control  input 
(Tschirha/t,  2003). 
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This  work  result^  in  a  legible  conirulltfr  aJ|(uriihin  assuming  perfeci  kAow)edge 
of  (be  lorgelsaiellice's  Slate  (posKion  and  velocicyL  The  flybrid  coolroUer  developed 
oebieves  rendezvous  (o  die  specified  relative  distance  and  velodry  in  590  mtnuies,  using 
4B.9in/6  AK.  The  TiOid  results  of  Tschlrtian'sconiroUeraiulysis  are  sumitiiirued  in 
Table  I. 


eo«i»^Type 


ImpiAfve  (CW) 

Conliiuous  <LQR) 
Hytnd 


~S^< 


;ca 

1  Tlm#te  wandMVDUi 

.C2; 

35  75 
3B3.11 


XS' 

3^4 

550 


*Mms  vnpuikKi^  ConteiiarOMa  twi  niMi  cnte<ia.  3  3  "m  ■  cjomw  eoproecn 


l(  is  ligniflcaol  lo  note  (be  Hybrid  cooiroUer  achieved  rendezvous  wt(h  considerable 
AK  savings  over  ibe  LQR  controller. 

Reconunendaiions  for  further  research  locluded  tovesiigaling  the  use  of  gam 
scheduling  as  part  of  an  LQR  controller,  and  (he  Incorpurauon  of  a  sequenddJ  filter.  Gam 
scheduling  was  suggested  in  order  lo  lower  conirol  u.sage  during  the  fnajon(y  ol'ihe 
rendezvous,  and  then  increase  tt  at  the  end  to  ci»nplese  ihe  reodezv'ous  without  (he 
complexity  of  a  hybrid  controller.  A  sequential  Qltec  was  recommended  lo  esiimaia  (he 
Slate  of  the  lorgei  saiclUle,  mcorporating  realistic  uncenamlies  in  using  sensor 
meosuremenis  (Tscbirban,  2003) 

The  LQR  controller  is  of  most  inieresi  to  this  researcher  given  the  gam-schedulmg 
recommeodattOD.  Therefore  (he  Dual  LQR  design  results  will  be  reviewed  here. 
Tschirbait  used  a  constant  State  Weighuog  Matrix,  Q  as; 
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I  0  0  0  0  0 
0  1  0  0  0  0 
0  0  10  0  0 
0  0  0  1  0  0 
0  0  0  0  1  0 
0  0  0  0  0  1 


iitvj  a  constant  Control  Wetgnimg  Matrix,  R  as; 


Rs 


Sen  0 
0  5^12 

0  0 


0 

0 

5«'12 


TtK  quadratic  cost  function; 


0) 


(2> 


(3> 


was  then  minimized  in  order  to  obtain  the  optinui  gain  matnx  to  oppiy  to  the  control 
thrust,  where  x  represents  the  system  state  (position  and  velocity!  snd  w  is  a  vector  of 
control  inputs. 

The  controller  decreased  the  relative  distance  between  the  microsatellite  and 
target  satellite  as  shown  below  in  Figure  S.  The  distances  in  the  figure  were  calculated  m 
the  relative  reference  frame  and  propagated  with  the  linear  equations  of  motion 
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Figure  5.  Relative  Di<iunce  during  LQR  Rendezvous  (atlee  Tschirhart,  2003) 

The  posiibon  ot'ihe  microsaielliie  relative  lo  ihe  target,  captured  m  ihe  dr,  r^66 
plane  IS  shown  in  Figure  6.  The  figure  nicely  illustrates  how  the  mjcrosatellite  initially 
begins  trailing  the  target  by  1000  km  in  the  same  orbital  plane,  and  then  drops  in  altitude 
to  increase  its  speed  (i.e  mean  motion).  Final  rendezvous  is  achieved  as  the 
fflicrosatelliie  arrives  within  1  m  and  1  cm/s  of  the  target  This  particular  controller 
configuration  led  to  the  final  LQR  rendezvous  results  included  in  Table  I  above. 
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Figure  6.  LQR  Rendezvous  in  the  Sr,r^S6  plane  (after  Tschirtiarr,  2003) 


Orbit  DeieriQloailniL/>lavlsoilon 

A  three-phase  tracking  system  architecture  concept  and  orbii  detemu  nation 
rouiinea  for  non-cooperative  rendezvous  were  developed  by  another  AFIT  master’s 
srudeni,  Bnan  Foster.  Of  particular  interest  to  this  research  is  the  on-orbit,  third  phase, 
orbit  determination  routine  to  estimate  the  target  satellite's  orbit.  A  Non-linear  Least 
Squares  orbit  determination  filter  was  implemented  to  accomplish  this  final  phase.  As 
expected,  the  filter  converged  to  a  solution  based  on  simulated  data 

The  orbit  determination  filter,  as  implemented,  was  found  to  perform  best  given  a 
large  number  of  observations  which  took  more  collection  time  and  thus  would  cause 
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significant  processing  delays  (Foster,  2003).  One  specific  simulation  run  found  the  filler 
was  able  to  reduce  the  estimate  error  from  an  iniiial  5.2  km  to  approximately  5  m  given 
100  data  points  (sensor  observations)  separated  by  60  seconds  each.  Less  data  soil 
allowed  ihe  filler  to  converge  on  an  estimate,  but  included  a  much  larger  error  compared 
10  die  truth  model. 

Foster  realized  that  in  a  rendezvous  mission,  time  to  collect  and  process  data  may 
noi  be  available  and  thus  control  maneuvers  may  have  to  be  based  on  less  accurate 
position  estimates.  The  develo^nent  of  a  Kalman-type  filter  to  allow  for  real-time 
processing  of  observation  data  for  die  oibit  determination  process  was  among  the 
recommendations  for  future  work. 

SuiDioary 

The  study  and  use  of  microsatellites  to  perform  a  variety  of  missions  is  currently 
underway,  (t  is  becoming  routine  to  not  only  consider  small  satellites  for  technology 
demonstrations,  but  operational  missions  as  well.  Industry  is  responding  to  demand  by 
producing  creative  solutions  with  applications  only  bound  by  human  imagination. 
Previous  AFIT  research  on  control  laws  and  orbit  determination  paved  the  way  for  an 
integrated  GN&C  analysis. 
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IK.  Mtftbodotogv 


Cbaplcr  Oven  lew 

The  purpose  of  (hi£  chopier  is  rwofold.  The  first  pan  describes  (be  sysieeis 
cnginecnnj;  appnuch  uken  for  this  feastbiltcy  analysis.  This  secuon  includes  o  lop-level 
systems  archiiecdire  fur  ihe  potential  system  bung  evaJuaied.  The  second  pan  deiails  the 
Qccesstiry'  lechnicaJ  theury  ret^uired  lo  solve  the  rendezvous  problem.  This  includes 
orbital  dyoiDiucs,  conirol,  orbii  deiermiDUdOD,  and  esiimaiioc  theory.  The  method  taken 
ui  A  top-down  systems  approach  with  the  majority  of  efibn  being  spent  on  the  dnving 
GN&C  algorithm  integrauon  Math  Works'  MATLAfi*  software  was  the  fool  used  for 
the  aJgunthm  development  and  evAluation. 

The  problem  statement  specified  that  the  microsatellne  will  begin  Approximately 
1000  km  bebmd  the  target  in  the  same  orbital  plane.  In  order  to  befior  scope  this  project, 
the  rendezvous  has  been  segineoied  into  phases  The  Dverardung  CONOPS,  Including 
the  three  phases,  is  in  Table  2  below. 


Table  2  Overarebutg  CONOPS 


fintf  urn 

! 

i 

Pwpgll 

OC-1 

1000  km 

1000  km 

Gmund 

Obtain  Intaal  Eatimsle 

OC-2 

1000  km 

S  km 

Gmund 

Ifxttal  RarMervoua 

OC-Z 

Sim 

1  m 

On-Otii 

Frtel  Rendezvous 

In  Overarching  CONOPS  Riase  l,OC-i,  ground-based  sensors  would  be  used  to 
generate  on  initial  target  state  estimate  Although  this  early  (riiase  would  likely  require 


signtficimi  global  mhastrucrure  to  do  well,  it  is  not  the  focus  of  this  research  (Foster, 
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2003)  Pha.ce  2  involves  closing  ihe  relative  dj stance  between  ihe  microsaieiliie  and 
largei  down  lo  S  km.  Tins  value  was  chosen  as  ii  represents  the  ex[tected  outside  range 
of  an  on-orbit  lidar  sensor.  The  control  law  work  ofTschirhan  led  this  researcher  to 
determine  this  second  pha.ce  is  quite  achievable  to  a  reasonable  error.  The  final 
rendezvous  phase,  OC-3,  is  what  is  studied  in  detail  in  this  work  The  control  and 
estimate  accuracy  required  to  achieve  rendezvous  to  within  I  m  is  certainly  the  most 
challenging  pan  of  the  problem. 

Systems  Engineering  View 

There  is  a  distinct  difference  between  traditional,  or  discipline-specific 
engineering,  and  systems  engineering.  According  to  Dennis  Buede,  a  well  respected 
evpen  in  the  systems  engineering  field.  Engineering  is  defined  as  a  “discipline  for 
translbrming  scientific  concepts  i  nto  cost-ellective  products  through  the  use  of  analysis 
and  Judgment”  (Buede,  2000).  This  often  applies  bent  to  hardware  component  or 
individual  software  item  development  Buede  further  defines  the  Engineering  of  a 
System  to  be  the  “engineering  discipline  that  develops,  matches,  and  trades  off 
requirements,  functions,  and  alternate  system  resources  to  achieve  a  cost-effective,  life¬ 
cycle-balanced  product  based  upon  the  needs  of  the  stakeholders”  (Buede,  2000). 
Systems  engineering,  at  a  very  basic  level,  is  the  effort  to  create  an  entire  integrated 
system,  not  jUst  a  bunch  of  components,  to  satisfy  the  need. 

Taking  a  systems  view  involves  the  up-5^3nt  planning  for  and  subsequent 
integration  of  the  traditional  engineering  products.  Many  standard  tools  are  becoming 
available  lo  the  systems  engineer  which  result  largely  in  non-matenal  products  essential 
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10  system  analysis  or  evaluation.  A  concepts  applicable  lo  the  analysis  of  a  system 
include  an  Operaiional  Concept,  External  Systems  Diagram,  and  a  Systems  Engineering 
Archiiecmre. 

An  Operational  Concept  often  includes  a  vision  for  wpat  the  system  is,  a 
statement  of  mission  requirements,  and  a  description  of  how  the  system  might  be  used. 
Figure  4  below  shows  three  primary  choices  considered  by  NASA  engineers  in 
determining  an  Operational  Concept  for  the  moon  landing  during  the  I960’s  (Brooks  et 
al,  1979;  Murry  and  Cox,  19891. 


Figure  7.  Alternative  Concepts  for  Apollo  Moon  Landing  (Murry  and  Cox,  19891 


This  illustration  demonstrates  how  several  potential  alternatives  may  exist  for 
solving  a  problem.  The  selection  of  the  most  desirable  concept(s)  is  the  Erst  step  in 


22 


eva)iintmg  the  feasibihry  of  the  system.  Cieeriy,  if  a  feasible  Operaiioiul  Concept  exists, 
then  It  is  possible  that  a  material  solution  can  follow  Examples  of  Operational  Concepts 
that  did  not  work  out  in  practice  include  those  for  prevKius  missile  defense  programs  such 
as  the  Strategic  Defense  Initiative,  Brilliant  Eyes,  and  Brilliant  Pebbles.  These  cases 
show  that  it  IS  not  sufficient  to  have  just  an  Operational  Concept  It  can  identify  flaws  in 
initial  thinking,  but  cannot  definitively  tell  you  the  system  will  work.  More  efton  is 
needed  for  that 

The  creation  of  an  External  Systems  Diagram  (ESDI  is  another  useful  tool  in  the 
design  and  evaluation  of  a  system.  It  is  a  meta-system  model  of  the  interaction  of  the 
system  with  other  external  systems  and  the  relevant  context  (Buede,  2000).  The 
recognized  value  of  an  BSD  is  in  clearly  defining  system  boundaries.  Although  these 
boundaries  have  many  useful  roles  for  the  systems  engineer,  for  this  project  they  simply 
help  put  the  system  in  context  to  aid  in  the  feasibility  asses^nent  An  ESD  can  be 
depicted,  IQ  Its  simplest  form,  as  in  Figure  ^  below.  The  system  itself  external  systems, 
and  the  context  can  all  be  clearly  differentiated  using  a  model  of  this  type 
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Context 


N 


External  Systems 


System  ''  ) 


are  impacted  by  ^'8y«tam^  \^J 

impacta.  but  r>ol  impacted  by.  **Syatam"" 


Figure  8.  Depiction  of  the  System,  EviBnul  Systems,  end  Context  (Wie^nge,  199S) 


A  Systems  Engineefing  Architecture  is  useful  for  creeling  (i.e.  concepruDliiing, 
designing  and  building)  complex,  unprecedented  systems  Architecting  is  known  to  be 
both  an  ait  and  a  science  in  both  the  traditional  home  building  and  space  system  domains. 
Architectures  are  not  jusi  useful  in  the  development  of  systems,  however.  An  emerging 
application  is  in  carrying  out  behavior  and  performance  analysis  and  to  evaluate  potential 
system  designs.  Specifically,  architectures  are  beginning  to  be  used  to  help  determine  if  a 
proposed  system  will  perform  the  desired  mission  in  the  desired  manner,  or  Operational 
Concept 

A  complete  systems  architecture  is  composed  of  the  three  views,  or  perspectives. 
Figure  9  below  illustrates  how  the  Operational,  Systems,  and  Technical  Standards  Views 
are  combined  to  fully  describe  the  system. 
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Figure  9  Arcluievture  Views  (DoDAF) 


All  archiieciure  deveiuped  fur  t  specific  purpose  und  only  to  ihe  poiAi  ibui  uwfu) 
results  are  obtained.  The  right  mix  orhigh-ieve)  and  deuil  views  must  be  ^ught  fur  an 
eOective^  eiUcieol  onitict  to  result 

The  Operational  View  (OV)  includes  the  tasks,  activities,  and  operauaoaJ 
elements.  It  generally  involves  both  graphic  and  teunjal  descriptions  to  convey  the 
concepts  and  intended  uses  of  the  sysrem  TheSysrons  View  describes  und  interrelates 
the  technologies,  systems  and  other  resources  necessary  to  support  the  requirements.  The 
Technical  Standards  View  contains  the  rules,  conventions  and  standards  governing 
system  imploneniation. 

The  intent  of  this  research  is  tu  develop  unly  the  minimum  set  uf  architecture 
products  necessary  to  make  a  top-level  evaluation  of  system  anbievabiliey.  This 
researcher  bas  developed  three  DV  products  tor  the  sysrem  evaluation:  High-level 
Operational  Concept  Graphic,  Operational  Concept  Narrative,  and  FuncuanaJ 
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Decomposinon.  Therefore,  ihe  Sysiems  ond  Technica)  Standards  views  have  no<  been 
corepieted  Alibougb  all  three  views  are  necessary  for  a  complete  systems  archilecnire, 
the  OV  IS  deenied  sulBeient  for  this  prejeci  evelimlioa. 

la  rormally  eveJuijimg  a  system,  the  development  uf  ibe  sysiem  arvhiieciure  ts  just 
one  step  m  the  process.  Dr  AJevonder  LevLs,  Chief  Scteniisi  of  the  Air  Force,  outlines  on 
evuluauon  approach  in  Figure  10  below.  Once  an  archtieccure  is  developed,  on 
executable  model  nuisi  iheo  be  consirucied  and  rua  lo  develop  analysis  results.  Only  a 
top-level  arcbiieccure  design  was  developed  for  this  project,  therefore  only  a  (^ualiiaiive 
ev  al  uaiion  of  the  sysiem  can  be  made. 


Figure  10.  Arcbiiccture  Evaluation  Approach  (afler  Levts,  2QQ3) 


Cuairul  Theory 

This  section  uuUines  the  oibiiul  dynaimcs  theory  applied  to  the  control  aspects  of 
the  rendezvous  problem.  Guidance,  ur  orbit  control,  is  delined  simply  by  Wertz  as 
“odjUiimg  the  orbit  to  meet  some  predeienmned  coadiuoos”  (Wertz,  1999).  The 
conditions  to  this  case  are  those  of  a  successful  reodezv'ous,  nonunally  witbm  )  m 
relative  distance  and  1  cm/s  relative  velocity  between  the  microsaielliteand  tui^eL 
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Before  the  details  S[>ecitic  to  this  rendezvous  [iroblem  ore  discussed,  an  interesting 
hisioricol  control  system  is  presented  Perfiaps  one  of  the  earii  est  control  systems  ever 
employed  was  by  Hero  of  Alevandna  in  ancient  times  The  device  for  opening  his 
temple  doors  is  shown  in  Figure  1 1  below. 


Figure  11.  Hero’s  Control  System  for  Opening  Temple  Doors 


The  system  input  was  lighting  the  alter  Tire.  Water  from  the  container  on  the  left 
was  dnven  to  the  bucket  on  the  right  by  the  expanding  hot  air  under  the  tire.  The  bucket 
descended  os  it  became  heavier,  thus  turning  the  door  spindles  and  opening  the  doors. 
EKimguishing  the  fire  had  the  opposite  efteci.  As  the  control  mechanism  was  not  known 
to  the  masses,  it  created  an  air  of  mystery,  demonstrating  the  power  of  the  Olympian 
gods. 
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In  order  to  apply  the  tiKory  (o  CM  rendezvous  problem,  ii  is  necessary  (o  first 
ouilme  relative  motion.  This  theory  describes  the  nucrosatelliie  and  target  positions  and 
velociiies  relative  to  a  circular  reference  frame.  Hill’s  coordinate  frame,  shown  in  Figure 
12,  can  aid  in  illustrating  this  concept.  The  origin  O  is  centered  in  the  Earth  and  fixed  in 
inertial  space.  V  is  the  origin  of  a  reference  ftame  that  is  centered  on  the  instantaneous 
location  of  a  point  moving  about  O  in  a  circular  orbit  with  mean  motion,  n .  The  unit 
vectors  in  the  circular  reference  frame  (RTZ)  are  m  the  radial,  in-track,  and  out 

of  plane  directions,  respectively,  and  r.  is  the  radius  of  the  circular  reference  orbit 
(Tragesser.  2003). 


Figure  12  Hill’s  iRTZ)  Coordinate  Frame 

A  satelliie  can  be  added  to  Figure  12,  to  illustrate  the  relative  position  from  the 
reference  orbit.  Figure  13  below  illustrates  this  relative  position  in  the  RTZ  coordinate 
h^me. 
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S:neLliie 


Figure  13.  Relative  Posiiion  m  RTZ  Coordinate  Frame 


In  thia  frame,  ihe  posiiion  of  the  feline  is; 

^  *  |(';  ^)C08^£>^,  +  ((f,  -K 

and  ihe  velociry  can  be  found  from' 


(41 


(5) 


where  die  superscn^iia  i  and  ••  correspond  to  die  inema)  and  circular  reference  frames, 
respectively,  and  ihe  mean  motion  of  the  circular  reference  frame  is' 


(61 


Through  a  fair  amouni  of  manipulailon,  the  relative  equations  of  motions  follow  as  in 
Equation  7. 
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£r  —  Irtr^Sd  —  n  *  (r  •«>  s  — {r,  —  2dt) 
r^Sd  *  2fl  df*  -  ff  *  r^S6  m  -n^  r^SB 


Solving  Equation  7  yields: 
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Equations  10  and  1 1  descnbe  ihe  relative  veloaty  and  position,  respectively,  of  the 
satellite  with  reference  to  a  circular  reference  orbit  For  the  theory  to  twld,  both  the 
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mtcrosaielliie  uid  target  musi  remain  suiTlaently  close  to  ihe  circular  reference  orbit 
(Wiesel,  1997). 


Tbe  specific  control  ler  scheme  chosen  uses  a  Linear  Quadratic  Regulator  (LQR) 
following  the  successful  resulis  of  Tschirhart's  work.  The  theory  required  for 
implementation  follows,  based  on  the  rebtive  reference  ftame  in  Figure  13  and  the 
relative  equations  of  motion  in  Equation  7. 

A  state  vector  comprising  the  relative  velocity  and  position  can  be  defined  as' 
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and  the  derivative  as' 
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(13) 


The  relative  equations  of  motion  can  be  placed  in  state  equation  form' 

(14) 

where  u  is  a  vector  of  control  inputs  Equations  7. 12  and  1 3  can  now  be  used  to  rewrite 
Equation  14  as: 
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A  Linear  Quadratic  Regulator  obuin&  the  optimal  gain  matrix  K  such  that  (he  stale- 
feedback  law; 

w»~AS  (I6i 

minimizes  the  quadratic  cost  function; 

^sf{x’0x-*-e/'J?n)i//  (I7> 

{I 

The  as»cidte^  Rh:cati  equation  is  solved  for  5 : 

SA&A'S-SBR  (I8i 

where  Q  is  the  State  Weighting  Mainx  and  B  is  (he  Control  Weighting  Matrix 

Higher  values  in  the  Q  matrix  speed  movement  toward  (he  desired  state,  and 
higher  values  in  the  A  matrix  reduce  control  usage  (Tschirhart,  2003).  The  values  of  Q 


and  R  have  been  selected  to  follow  the  forms  of  Equations  1 9  and  20  below  The  values 
of  If  were  he  set  to  1,  while  r  was  allowed  to  vary  during  (he  rendezvous  process  as  gain 
scheduling  is  implemented. 
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MATLAB’s  LQR  tUnciion  is  used  lo  calculate  ihe  opiimal  gain  tnairiv  as; 


K  sj?  '5*5 


(I9i 


(20) 


(2\) 


TlK  conirol  mpui  of  Equation  16  must  now  be  modified  lo  account  for  ilK  fact  the 
microsatellite  is  chasmg  IM  target  rather  than  the  reference.  This  control  should  be  based 
on  the  difference  between  the  microsaiellite's  state  vector  and  the  target's  state  vector' 
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(22) 


The  LQR  routine  developed  by  Tschirhait,  outlined  m  Figure  4,  was  used  to  begin 
this  study.  The  final  hybrid  controller  solution  was  not  examined  in  favor  of 
implementing  gam  scheduling  in  the  LQR  algorithm.  Once  better  understood,  this  gam 
scheduling  LQR  conuxiller  was  coupled  with  different  estimation  filters,  attempting  to 
construct  an  integrated  solution 
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Orbit  Deicrmlnailon  Theory 

In  using  only  the  above  orbsial  dynamics  and  control  theory  to  solve  the 
rendezvous  problem,  perfect  knowledge  of  both  the  nucrosatelhie  and  the  target  must  be 
assumed.  One  can  expect  the  microsaiellite  maintains  ^trly  good  knowledge  of  Its  own 
stale,  by  using  GPS  for  example.  The  position  and  velocity  of  the  target,  however,  must 
be  estimated  in  some  manner.  As  in  an  electrical  filter  which  extracts  the  desired  signal 
^om  the  undesired,  an  estimation  algorithm  which  extracts  the  system  state  from 
observations  with  errors  is  called  a  filter  ( Wiesel,  2003).  An  observer,  or  filter,  must  be 
designed  to  estimate  the  plant  states  that  are  not  directly  observed. 

There  are  various  methods  of  reconstructing  the  states  from  the  measured 
outputs  of  a  dynamical  system.  Estimation  tiller  types  relevant  to  this  research  include: 
linear,  nonlinear,  batch,  and  sequential.  A  linear  estimator  assumes  the  data  is  linearly 
related  to  the  system  state  at  the  time  taken.  This  can  greatly  simplify  the  problem,  but  is 
not  applicable  in  many  cases  In  a  nonlinear  estimator,  the  observed  quanuties  are 
allowed  to  be  related  to  the  system  state  by  a  very  nonlinear  set  of  relations  ( Wiesel, 
2003)  Another  way  to  classtly  a  filler  is  via  how  the  data  are  processed.  A  batch 
algorithm  assumes  all  data  are  available  before  the  estimation  process  begins,  and  is  all 
processed  in  one  large  batch.  A  sequential  filter  is  continuously  taking  in  new  data  and 
producing  an  improved  estimate.  A  batch  algorithm  can  be  made  more  sequential  by 
observing  and  processing  smaller  batches  of  data. 

An  estimate  of  die  target  state,  i ,  will  be  represented  as  x .  A  trajectory  to 
linearize  the  dynamics  about  must  also  be  chosen.  As  the  true  system  state  x  is 
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imobuinjble  and  an  esiitnaie  i  doe&  not  yet  exist,  a  reference  trajectory  ts  used  The 
reference  trajectory,  ,  is  a  trajectory  one  expects  will  be  close  to  the  estimate.  The 
goal  IS  to  find  correciions  to  the  reference  trajectory  turning  it  into  a  reasonably  good 
estimate,  x .  The  reference  trajectory  usually  comes  from  an  initiaJ  orbit  detemunation 
method  and  is  then  updated  based  on  observation  data  ( Wiesel,  2003). 

The  control  law  in  Equation  16  will  be  modified  asi 

5  »  -AS  (23) 

For  a  linear  filter,  a  reasonable  way  to  estimate  i  is  by  duplicating  the  actual  state 
dynamics  in  propagating  i  (Cobb,  2003).  The  equations  of  motion  for  the  target  as 
shown  in  E4^uation  14  will  then  be  described  as; 

(24) 

Since  the  true  target  state,  x  cannot  be  measured,  a  correction  term  must  be  added 
to  the  dynamics  equation  for  the  observer.  Equation  24  can  be  modified  to  include  a 
correction  term  proportional  to  the  difference  between  the  measured  and  estimated  output 
(Cobb,  2003) 

f9A!*Su*l{y^Cy)  (25) 

where  L  is  the  estimator  gain  matrix,  y  represents  the  mode)  for  the  measurement,  and 
C  models  the  observation  geometry.  The  y^Cl  term  represents  the  filter  residuals. 

This  residual  term  can  also  be  described  as  the  difterence  berween  measured  and 
predicted  observation  values; 


35 


(26i 


Shown  in  the  form  ofE^^uation  26,  n  is  ovideni  the  goal  of  ihe  filter  is  lo  minimue  the 
resHjua),  or  correction  term. 

MATLAB  contains  a  number  of  built-in  linear  filters.  A  sim^ile  one  lo  use  is 
LQE.  The  Linear  Quadratic  Estimator,  or  LQE,  follows  a  stationary  Kalman  estimator 
design  for  continuous-time  systems,  (t  returns  the  observer  gain  matrix  L  to  include  in 
the  equation  of  motion  given  in  Equation  25.  A  block  diagram  showing  how  a  Kalman 
filter  can  used  to  form  a  Kalman  estimator  is  shown  in  Figure  14. 


Figure  14  Kalman  Estimator  (MATLAB) 


The  MATLAB  function  LSIM  can  then  be  used  to  simulate  the  model  response, 
obtaining  the  state  estimate  x  Specifically,  the  command  Isim  (sys,  u,  t)  produces  a  plot 
of  the  time  response  of  the  LT(  model  sys  to  the  input  time  history  t,  u. 

As  stated  above,  a  simple  linear  filter  is  not  sufficient  for  many  control  problems. 
In  such  cases,  a  more  complex  nonlinear  estimator  may  be  required.  'To  handle 
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problems  which  are  useful  m  ihe  real  world  we  musi  abandon  the  linear  case  and  work 
with  nonlinear  system  dynamics  and  ihe  nonlinear  obtervaiion  geomeiry”  (Wiesel,  2003) 
Wiesel  details  a  very  good  algonihm  for  a  Non-linear  Least  Squares  esiimaior  m  his  book 
Modern  O-iiif  Delennmaunn.  This  routine  amounts  to  calculating  the  state  from  the 
observations.  The  specific  steps  given  by  Wiesel  are  as  follows: 

I  Propagate  the  state  vector  to  the  observation  time  i,  and  obtain  the  state 
transition  matrix  ) 

{ ds  ^  dx  cii  dx  dx 
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and  x,y,2  refer  to  the  three  components  of  the  position  vector.  The  state 
transition  matrix  comes  from  linear  dynamical  systems,  where 
®  propagates  the  actual  state  as  a  function  of  time.  It  is  the  gradient  of  the 
solution  with  respect  to  the  initial  conditions 
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Obtain  the  residual  vector  r,  s  r,  -  G(z) ,  where  r,  is  the  measured 
observation  vector  and  G(i)  is  the  predicted  data  vector  of  the  current 
state  vector  x .  Calculate  the  observation  mode!  /f,  for  this  particular 


data  point,  where  /f. 


Then  calculate  the  observation  matrix 
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Add  new  terms  to  the  running  sums  of  the  matrix  ‘t,  and  the 

vector 
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where  Q  is  ihe  coral  insinimeni  covariance  mainx 


Whan  ail  dala  has  bean  pm<.‘&ised: 

4  Calculate  Ihe  covanance  of  the  correciun  'TJ  '  and  (he 

stare  cotrecrion  vector  ar  epoch  SrU, ) » /  «  ^  7]  y,  'r,  Note  tliac  the 
matri  x  T,Q''T,  tnusi  be  invertible  for  a  new  esiimate  of  the  reference 
trajectory  (o  exist  This  is  known  as  (he  observabiliiy  condition. 

5  Correct  the  reference  ira^crory  )  •  •‘rrr.i 

new  estimate  of  the  reference  trajectory 

6  Deremune  if  (he  process  has  converged.  If  not  begin  again  at  step  1.  If 

80,  18  an  estimate  with  covariance  P, . 

7  Check  to  ensure  there  are  no  unbelievably  large,  greater  than  3^ , 
residuals.  If  so,  reject  (he  observation  in  step  3. 

Typical  observation  measurements  are  of  the  form'  range;  range  and  range  rate; 
range,  azimuth  and  elevation.  The  measurements  chosen  for  this  research  include  range, 
azimuth  and  elevation.  The  observations  r,  can  (hen  be  described  as' 


1* 

uzifnuth 

p 

a 

477777 

tan*'  {y/x) 

election 

P. 

tan  'Ir/Vr'  -y 

As  the  predicted  observation  value  G(i)  takes  the  same  form  as  ? ,  (he  observation 
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components  are  given  by: 
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(291 
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and  a)  I  other  components  equal  to  2efo.  H,  is  found  by  ukmg  the  partial  denvaiives  of 
the  G  vector  with  respect  to  the  state,  evaluated  on  the  reference  trajectory. 

In  non-linear  Hlter  design,  one  also  needs  to  decide  how  sequential  to  make  it 
The  work  of  Foster  was  based  on  a  strictly  batch  method,  assuming  all  data  are  available 
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and  prucHsed  at  once.  This  gencraOy  fine  for  many  applicattona  For  the  reodezvou^ 
problem  at  band,  however,  ii  would  bkely  take  loo  lon^  to  accuenulaie  al)  the  deaired 
obaervauon  data  and  proceaa  before  moklo^  an  updaie  (o  (be  target  csuenaie.  Wailing  for 
thi^  baicb  prooeas,  the  unoeruiniy  tn  the  estimate  oonusuvs  to  grow  os  (he  lorgei  and 
rmeroaacelliie  orbii  (he  earth.  Care  mu be  taken,  buwever.  to  nui  dvaign  a  "ily- 
folluwer"  which  may  atiempi  to  produce  a  new  cstmuie  given  only  a  few  obaervadon 
The  above  diacua&ion  relates  lo  waicluog  (he  covariance  of  (be  estimate,  which 
grows  during  propagation  It  also  lends  to  gel  smaller  given  new  data.  This  effect  is 
shown  in  Figure  IS.  where  the  slate  covariance  grows  between  updaies  and  drops  when 
updated 


Figure  15.  Covariance  Behavior  with  Time  (Wiesel,  2003) 


Ifthe  niter  behavior  Lsusihe  lopline  *To  Ignuronce."  it  foils  as  the  knowledge  ofibe 
curreni  system  stale  becomes  mure  and  mure  unoenain.  If  behaving  as  “To  PerfecKon.’' 
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ihe  filler  believes  ii  has  achieved  pe^ecQon  uid  will  cease  lo  update  the  esiimaie  even 
given  new  observation  data.  This  condition  is  olfen  referreij  to  as  smugness  as  the  filter 
wi  II  not  react  to  additional  input.  Neither  of  the  above  two  conditions  is  desirable  and 
avoiding  them  is  an  art  of  filter  design. 

SuiDioary 

This  chapter  outlined  the  Systems  Engineering  approach  taken  to  evaluate  the 
non-cooperative  rendezvous.  An  Operational  Concept,  External  Systems  Diagram,  and 
Architecture  were  described  as  tools  to  assess  top-level  system  feasibility  A  Linear 
Quadratic  Regulator  was  discussed  in  context  of  orbit  control  laws.  Both  linear  and  non¬ 
linear  filter  theory  was  given  to  estimate  the  target  state  Wiesel's  Non-Linear  Least 
Squares  algorithm  was  detailed  as  a  speufic  filter  routine.  Finally,  a  few  subtle  filter 
design  considerations  were  discussed. 
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(V.  Analysis  and  Results 


Chapter  Overview 

This  chapter  contains  the  developed  systems  engmeeeing  products,  techmcai 
ON&C  algorithm  analysis,  and  subsystem  integration  results.  The  systems  architecture 
and  associated  products  developed  are  a  necessary,  but  insufTIcient,  step  in  the  feasibility 
evaluation  Assessment  of  the  select  products  indicates  top-level  system  feasibility  while 
underscoring  technical  and  integration  complexity. 

Beginning  the  technical  study,  previous  control  law  development  was  extended, 
by  gain-scheduling,  to  show  positive  trade  space  between  lime-to-rendezvous  and  fuel 
usage.  Given  this  positive  result,  attention  was  focused  on  the  orbit  determination  filter. 
The  use  of  a  linear  estimator  is  shown  to  be  inappropriate,  while  a  nonlinear  esiimaior 
requires  advanced  implementauon  for  the  application,  (niegraiion  of  tailored  controller 
and  estimator  components  proved  to  be  beyond  the  limits  of  text  book  algorithms.  It 
would  be  a  non-trivial  task  to  improve  these  algorithms  to  account  for  the  necessary 
complex  orbital  dynamics  involved. 

Extending  this  result  leads  to  a  low  probability  of  designing,  building,  and 
operating  a  microsaielliie  to  rendezvous  with  non-cooperaiive  targets,  using  established 
ON&C  software  routines,  in  the  very  near  term.  A  systems  or  technical  view  is 
insufficient  to  show  feasibility  by  itself  The  details  that  follow  show  how  one  view  leads 
to  possible  aiiainmeni,  while  the  other  points  to  serious  challenges. 
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Sykiem^i  Enuinvertng  Froai-ead  Ruulti 
Opvraiioool  Concept 

There  &re  a  oumber  of  waya  to  which  a  aaielliie  can  perform  a  non-cooperative 
rendezvoui  with  a  lar^t  aaieDiie.  The  chase  saielhte  c'an  be  Jirecily  launched  to 
rendezvoui  or  h  ciui  perform  orbit  traiufcrs  from  a  aimslar  orbK  as  the  target.  In  the  Inner 
ca&e,  trades  are  avn^lnbie  between  on-orbii  and  ground  sensor/proces^uig  ocuvities. 

Three  aliemsuve  Operational  Concepts  ore  depicted  m  Figures  16-18  below. 


Figure  16  Launch  to  RenJezvoiLs  Ops  Concept 


The  launch  to  rendeTVOus  concept  depicted  in  Figure  Id  above,  involves  critically 
deiermiDing  the  launch  irnuogoftbe  chase  saielliie.  The  most  efficient  way  to 
rendezvous  in  this  miinner  is  to  launch  when  ihc  target  orbit  passes  directly  over  the 
Launch  site.  Any  deviation.^  from  the  perfeci  Launch  tunc  cun  result  is  very  costly  orbit- 
pLooe  changes.  The  costs  are  not  only  in  required  weight  for  propulsion  lUel,  but  also  to 
design  complexity.  This  concept  is  relatively  compbcaied  and  difficult  to  execute,  in 


43 


order  lo  acuompli&h  the  rendezvous  dstecily,  wiihoui  any  on-orbit  maneuvers,  extreme 
Louoch  precision  h  required  This  further  requires  very  precise  knowledge  of  (he  taigei 
Siam  Given  the  above  prereqmsries,  (he  direct  launch  lo  rendezvous  operauooaJ  concept 
may  noi  represent  a  luw  lech  solution  very  well 


Figure  1 7  Autonomous  Rendezvous  Ops  Concept 


Fur  the  autonomous  rendezvous,  Figure  17.  the  chase  satellite  would  perform  the 
neeesstiry  orbit  determination  and  eonirol  activities  once  given  an  initiuJ  target  state 
estimate  from  the  ground  segment  Radar  and/or  optical  ^und  satellite  tracking  stations 
would  be  required  to  perform  initial  orbit  deierminaiion,  bui  then  yield  to  on-orb  it  sensors 
once  available.  Nominalty,  only  the  final  5  km  of  rendezvous  would  be  performed 
completely  autonomously  by  the  chase  satellite  This  is  primarily  limited  by  sensor 
performance  charactensiics.  The  autonomous  rendezvuus  concept  stilJ  requires  launch  to 
a  similar  orbit  as  the  target.  Just  not  requineg  (he  exactness  of  (he  direct  launch  concept 
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As  long  os  iheuhoie  vehicle  is  piJceU  tnio  lheu(neod^ilaJ  plane,  tnu ling  the  largtiby 
some  reasonable  amoimi,  (he  auiooomaijs  rendezsous  coocepi  retains  gre^i  poteolbl 


Figure  IB.  Ground-Assisi  Rendenous  Ops  Concept 

A  coocepi  using  ground  sensors  and  algoriihms  for  (he  eiuire  rendezvous  is  shown 
In  Figure  18.  Whereas  (he  previous  concept  only  used  ground  resources  for  (he  iniiiaJ 
orbit  deienruoaKon,  this  scheme  relies  on  ground  input  Tor  ihe  enure  rendezvous 
se<;uence  This  involves  eviessive  development  and  infrastructure  on  ihe  ground,  bu! 
relatively  biiie  on  die  chase  saieliiie.  The  saielliie  would  receive  specdfic  GN&C  and 
propulsion  commands  to  execute  each  rendezvous  maneuver  This  aliemauve  may 
alleviate  some  of  (he  spacecraft  development  challenges.  The  comprehensive  ground 
infrastructure  required,  however,  may  also  be  beyond  a  low  teeb  approach. 

Recalling  from  die  problem  suiemeni,  the  microsatellite  is  set  lu  begin  in  (be 
same  orbitaJ  plane,  approximately  1000  kna  behind  die  target  satelJiic.  The  microsaiellut! 
(hen  performs  rendezvous  maneuvers  to  approac*h  the  (orgeL  Therefore  the  autonomous 
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and  ground-asstsi  approaches  seem  to  be  more  appropnuie  Wan  We  launch  lo  rendezvous 
concept.  The  lanet  couJd  be  reviiiied,  however,  if  other  aJiemaiives  faded  (o  produce  an 
acceptable  soluuon.  The  above  work  establisbes  that  there  is  adequate  trade  space  to 
more  fully  define  a  final  OperutiondJ  Concept  and  thus  a  maierioj  solution  may  exist 

Evieroal  Systems  Diagram 

The  first  step  token  in  developing  an  ExteroaJ  Systems  Diagram  was  to  determine 
the  first-order  functions  of  the  sysiefo.  DeCmng  the  top-level  function  of  the  sysurmto 
be  ''Rendezvous  with  Target,”  three  additional  sub-functions  were  determuiciL  The  hi^- 
level  functioaal  decomposition,  in  Figure  1 9  below,  shows  acquiring  the  target  state,  and 
determining  and  executing  AF  maneuvers  as  sub-funciions  This  simple  decomposition 
contains  Functions  allocated  principally  to  propulsion  and  GN&C  subsystems 


Figure  19.  High-Level  Funciional  Decomposition 

The  standard  process  modeling  technique,  IDEFO  (Integrated  Definition  for 
Function  Modeling),  was  used  to  develop  the  External  Systems  Diagram  (ESD).  The 
decomposiiioD  above  was  used  as  a  stoning  point  to  begin  the  process.  Each  of  the  three 
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sub-iunctions  was  iJeniified  as  an  aciivliy  in  ihe  ESD.  Input,  uuipub  and  ounirols 
(triggers)  were  then  esiabbshed  for  eaiA  aclivily  fuociioo  The  resisting  Jiagnun  in 
Figure  20  definitively  displays  the  iraosformaiion  of  inputs  into  ouipuii  by  the  sysien. 


Figure  20.  R^dezvous  Euemal  Systems  Diagram 


The  system  boundary  is  not  yet  enpliuitly  clear  from  Figure  20.  however. 
TherEfore,  a  second  system  ESD  was  developed  following  the  vonsirtici  of  Figure  8.  The 
system  as  reprsenicd  ic  Figure  21  is  belter  dilTereoitaied  from  the  enviroiutiefii.  both 
external  sysiems  and  context  In  this  diagram,  ground-based  sensors  and  compuuog 
resources,  support  systems  and  ihe  target  satellite  are  depicted  as  exxemaJ  systems  which 
inieraci  with  the  microsdielliie.  The  context,  space  environment,  acts  on  the 
microsaielliie  system,  bm  i.s  clearly  noi  impacted  by  the  system. 
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Corttavl 


Even  though  boih  Exterrul  System  Diagrams  developed  are  relatively  simple,  they  do 
help  scope  the  system  and  are  a  necessary  prerequisite  to  system  architecture  creation. 

Systems  Engineering  Architeciure 

A  few  top-level  systems  architecture  products  have  been  developed  to  illustrate 
one  possible  mjcrosatellite  architecture  to  achieve  a  non-cooperative  rendezvous.  The 
products  include  an  Operational  Concept  Ciraphic,  accompanying  Narrative,  and  a 
detailed  Functional  Decomposition.  Taken  together,  they  document  a  rough  Operational 
View  of  the  architecture  As  stated  in  the  methodology,  only  a  qualitative  assessment  can 
be  made  based  on  these  products.  The  necessary  executable  models  to  base  a  quantitative 
assessment  on  are  well  beyond  the  scope  of  this  research. 
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Opvrutionul  Coocvpt  Graphic 

The  High-Lev  e)  Operauooal  Concept  graphic  for  (he  oucrosaiollKe  system  is 
shown  in  Figure  22.  The  Operational  Concept  Norrauve  that  follows  desenbes  in  more 
detoii  (lie  inieractions  between  the  entities  portrayed  in  the  gruphic. 


OperaiioDoi  Concept  Narrative 

The  researcher  considered  a  key  thread  to  describe  the  OperutionaJ  Concept  •  the 
se<;uence  of  aciiviitcs  that  mke  piace  when  a  user  authorizes  the  Microsatellite  System  to 
rendezvous  for  (he  purpose  of  disabling  the  T argei  sateUitc.  The  architecture  is.  at  this 
poent,  proceedmg  forward  assuming  the  Autonomous  Reodezv'ous  Ops  Concept 
olienuiive  depicted  m  Figure  i  7 
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In  order  to  further  bound  die  system,  die  key  externa)  systems  and  interna)  system 
components  need  to  be  identified.  The  inteenal  system  components  of  the  Microsatelhte 
System  include  the  items  widim  the  boxed  part  of  the  graphic.  They  include  the  GN&C, 
Propulsion,  Associated  Subsystems,  and  System  Operators.  GN&C  and  Propulsion 
subsystems  are  called  out  separately  from  other  “Associated  Subsystems”  due  to  their 
strong  relationship  to  the  activities  in  die  Rendezvous  ESD  of  Figure  20. 

System  Operators  are  also  drawn  inside  die  box  due  to  dependency  on  them  for 
non-autonomous  operations.  It  ts  tempting  to  exclude  die  personnel  from  the  system  and 
only  evaluate  the  hardware.  The  researcher  has  elected  to  create  the  OV  architecture 
products  including  operators  to  provide  a  better  system  evaluation. 

Key  external  systems  include  the  Target,  Command  Authorities,  Ground  and 
Support  Systems.  The  graphic  depicts  how  these  systems  internee  with  die 
Microsatelliie  System  to  include  nominal  activities. 

The  rendezvous  process  can  be  described  as  follows:  Command  Audioniles 
determine  a  Target  i^uires  disabling.  This  rei^uiremeni  is  sent  to  the  Microsaiellite 
System  Operators  via  a  Space  Tasking  Order  (STO).  The  System  Operators  send  a  target 
acquire  command  to  the  Microsatelliie  Associate  Subsystems.  The  Microsatelliie  then 
acquires  the  Target  stale  estimate,  initially  via  Ground  Systems.  The  Microsatelliie  then 
begins  die  rendezvous  maneuver,  or  AK  calculations  Once  determined,  the  thrust 
requirements  are  sent  to  the  Propulsion  system  to  execute.  The  target  state  acquisition 
through  propulsion  thrust  loop  continues  based  on  Ground  System  tracking  until  the 
Microsatelliie  is  within  range  of  its  on-board  tracking  sensor. 
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Once  wiihm  range,  nominally  S  km,  (he  Micro&aielliie  swiccMs  lo  autonomous 
mode  for  closed-loop  tracking  and  control.  The  target  acquisiiion/propulsion  thrust 
sequence  continues  as  before,  but  now  based  on  space-based  sensor  measuremenis,  uniil 
rendezvous  is  achieved.  When  rendezvous  is  achieved,  Associate  Subsystems  are 
engaged  lo  disable  the  identiHed  target.  The  disabling  activity  is  outside  the  scope  of  (his 
architecture  and  thus  will  not  be  detailed.  Finally,  System  Operators  generate  a  mission 
report  for  Command  Authontiea. 

Functional  Decomposition 

Given  the  Eviemal  Systems  Diagrams  and  Operational  Concept  presented  above, 
ic  IS  now  possible  to  further  detail  the  initial,  high-level,  functional  decomposition  shown 
in  Figure  19.  Whereas,  the  high-level  decomposition  facilitated  ESD  development,  a 
more  detailed  decomposition  is  required  lo  better  evaluate  the  potential  architecrure.  The 
detailed  functional  decomposition  is  shown  in  Figure  23. 
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Figure  23.  DctalkJ  FunciiooaJ  Dccompo&inun 


OD)y  TWO  of  ibe  three  ini  tial  sub-fiinciion&,  Aci^uire  Target  Stale  aad  Deienmoe 
dV  Requirements  were  f^ber  decomposed.  The  ibird  sub-funciton,  Perform  AV 
Maneuvers,  did  not  require  such  due  lo  the  retauveompliciiyorihls  function  compared 
to  die  other  two.  Each  of  the  lowesi-iier  iuncuon  could  again  be  decomposed,  but  is  not 
necessary  for  the  purpose  of  Jus  architecture  evaluation 

Analyzing  ibe  fuDciiooal  decomposnion  of  Figure  23,  two  conclusions  can  be 
drawn  First,  tbe  top-level  function.  Rendezvous  with  Target,  quiie  easily  decomposes 
into  logical  sub-runciions  of  necessary  depth  and  breadth  to  create  a  more  detailed 
architecture.  This  lends  suppoti  to  system  feasiblUty  ai  a  very  broad  level.  Second,  by 
onolyaDg  the  lowest-level  sub-functions,  it  becomes  apparent  ibat  a  collection  of 
complex,  iDiegraied  activiiies  must  take  place  to  perform  The  rendezvous  mission.  The 
two  functions  dMomposed  to  the  lowest  level.  Acquire  Target  Stale  and  Determine  AV 
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Requirements,  are  [ihysically  allocated  lo  the  ON&C  subsystem.  Tbe  following  sections 
provide  the  results  of  detailed  CiN&C  algorithm  analysis. 


Cootroller  GtlO'SchedulIng  Results 

The  first  technical  area  investigated  was  improving  the  controller  developed  by 
Tschirhan.  Although  the  impulsive/continuous  bum  hybrid  model  achieved  a  relatively 
efTlcient  rendezvous,  implementing  a  gain-scbeduling  LQR-only  controller  offered 
potential  improvement  with  less  complevity.  Thus,  the  State  and  Conucil  Weighting 
matrices  of  Equations  19  and  20.  respectively,  were  examined 

To  achieve  more  efficient  control,  one  desires  less  control  (larger  R  or  smaller 
Q)  initially,  and  more  control  (larger  Q  or  smaller  )  as  the  relative  distance  is  reduced. 
This  allows  larger  df'  maneuvers  to  achieve  final  rendezvous,  while  avoiding  such 
costly  maneuvers  early  in  the  process.  Only  one  weighting  matrix  needs  to  change  to 
attain  the  desired  result.  The  Q  matrix  was  (ef(  constant  while  the  R  matrix  was 
programmed  to  decrease  over  time.  The  Control  Weighting  Matrix  was  set  to  vary  as' 
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and  lilsi  _  miH'  is  ih«  reloiive  diuance  beivp'ren  ihe  two  sdieUiies  This  pan  of  the  siuJy 
involved  varying  buUi  ihe  fuial  magnitude  of  r  ( ).  unJ  R  _  Iwivr . 


Reasonjbic  imde  space  discovered  between  die  AT  rec|iiired  and  die  Time  to 


Rendezvous.  The  resiilLs  ore  suinmorLeed  below  in  Table  3. 


Table  3  Gain-Sehedulint'  Controller  Results 
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The  Ttoie  to  Rendezvous  and  die  AK  required  have,  ui  general,  an  inverse  relauooslup. 
This  IS  due  to  the  rendezvous  tune  primarily  beuig  driven  by  ibe  Control  Weighting 
Matrix,  R.  Asa  larger  R  causes  rendezvous  to  be  achieved  quicker,  more  fuel  i» 
consumed.  This  correlation  is  iDusirated  in  Figure  24. 
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Figure  24.  Gam-Scheduling  Trade  Rjsuks 


The  Gain-Scliedubn^  trade  resulted  in  a  confmDer  that  can  be  both  more  fuel 
elUdeni  and  obtains  rendezvous  quicker  ibaa  the  hybrid  controller  developed  by 
Tschirhart,  while  tnainiaining  flejubiluy.  Case  l-B,  for  example,  achieved  rendezvoin;  in 
572  minutes,  using  34  99  m!t  AK .  This  case  exhibited  Fuel  perrormance  slightly  better 
than  the  Impulsive  thrust  and  Hybrid  Controller,  and  an  order  of  magnitude  better  than 
the  Conuauous  Non-Oain-Seheduied  LQR  model  shown  in  Table  1.  Tbe  Tine  to 
Rendezvous  was  sbgbtly  better  than  the  >lybnd  Controller,  while  takug  a  few  additional 
earth  ort»its  over  the  Non-Ciain-Scheduled  LQR  ControQer.  The  related  performance 
results  are  given  below  in  Figures  25  and  26. 
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Figure  2S.  Rebtive  Dtsience  during  LQR.  Rendezvous  •  Case  l-B 
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m  •MO)  «e  ««  401  •»!  D  too 

Figure  26  LQR  Rendezvous  in  ihe  6t\r,S0  pbne  Case  l-B 

Although  the  results  in  Table  3  are  only  a  subsec,  they  show  adei^uaie  trade  space 
exists  between  the  Time  to  Rendezvous  and  AV  rei^uiremenis  The  Baseline  Case 
represents  (he  Non-Oain-Scheduled  LQR  solution  developed  by  Tschirhart.  The 
subsequent  cases  all  used  some  level  of  gain  scheduling,  as  listed  in  Table  3,  to  alter  the 
dependent  variables.  Higher  values  of  and  R_  factor  chan  are  listed  were  found  to 
either  not  converge  or  ac  least  take  too  long,  making  (hem  essentially  not  constructive. 
The  main  MATLAB  routine  is  included  m  Appendix  A 
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Cooiroller/Esiliatior  Integrotion 


Linear  Filter 

A  simple  linear  estimaior  was  attempted  to  iniegraie  with  ihe  Oain-Schedulmg 
Controller  developed  above,  (n  order  to  first  acquire  an  appreciation  for  the  MATLAB 
LQE  and  LSIM  tools,  a  very  straightforward  estimator  was  developed.  The  code  is  in 
Appendix  B.  For  this  model,  the  equations  of  motion  in  Equation  1 5,  with  tt>e  control 
vector  set  to  zero,  and  an  Observation  Oeometry  Matrix,  C  containing  range,  azimuth, 
and  elevation  were  used  to  obtain  the  estimator  gain  matrix,  L  by  LQE  command  LSIM 
was  then  used  to  calculate  the  target  and  target  estimate  states  over  time  This  is  an  open- 
loop  simulation,  of  the  target  only,  to  characterize  the  perlbrmance  of  a  simple  linear 
estimator. 

The  first  simulation  began  with  the  target  state  estimate  and  target  state  equal. 
That  IS,  target  knowledge  error  was  set  to  zero.  The  iniiiai  target  state  vector,  in  the  form 
of  Equation  1 2  was  given  by; 
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The  simulation  was  run  for  10,000  seconds.  The  resulting  target  position  and  velocity 
errors  are  plotted  in  Figures  27  and  28.  Note  the  state  estimates  appear  to  track  the  truth 
in  ail  cases  throughout  the  simulation.  Data  taken  at  the  end  of  the  simulation  show  the 
total  target  position  error,  6.16e-0O4  m,  and  the  total  velocity  error,  1.76e-006  m.  confirm 


the  vi<;ua(  ;:anc1iL'<H)n.  The«e  valuer  are  reasonably  small  as  expected,  given  ibe  perfect 
imiial  guess  for  die  stale,  lo  tai:i,  ibe  values  are  indisuoguisbablc  in  die  dgures  below. 


Figure  27.  LQE/'LSIM  Target  Posiiton  Error  with  Perfect  liulia)  Estuuate 
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Figure  28.  LQE/LSIM  Taryet  Velocity  E«of  with  Perfect  (ntibi  Estimaic 


G)ven  ihe  e^tinuiur  wih  not  hkely  begin  with  a  perleci  initial  gues^,  a  I  km  target 
position  error  was  loiroduced  The  1  km  value  was  assessed  as  the  best  an  initial  orbit 
deienninaiiOD  method  could  produce  as  a  starling  point  for  any  on-orbit  estimaiion 
(Foster.  2(X)3)  Using  the  same  procedure  as  above,  ihc  target  knowledge  error  was  set  to 
I  km.  split  evenly  between  the  e,  and  components  Tbe  li  Iter  performance  fortius 
caw  IS  given  in  Figures  29  and  30 
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Figurt  29,  LQE/LSIM  Target  PoiHion  Error  wiift  I  km  tnnial  Error 
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Figure  30.  LQE/LSfM  Targei  Velocity  Error  with  I  km  Initial  Error 

Although  the  suie  eBiinuie»  for  the  I  km  imual  error  caw  also  appear  to  track  the 
iruib  &Qmev>'hai  well,  a  oloser  inspeciion  is  required  Li  this  case  tl  is  possible  lo 
distiDguisb,  graphically,  between  the  truth  and  estimate.  la  Figure  31  bebw.  the  theta 
component  has  been  removed  to  better  illustrate  the  r  and  :  components  of  the  posiiloa 
truth  OJtd  estimate. 
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Figure  3 1 .  LQE/LSIM  Target  Pusaiijn  Error  »iUi  I  lun  IniitaJ  Error,  r  uiJ  ;  only 

Dciu  Liken  qi  ihe  erhJ  of  ihe  simulation  show  the  lou)  targ«it  po^jiion  error,  m,  and 

the  laai  vdociry  error,  1.40  m,  are  signilicimL  Therefore,  the  filter  does  not  exhibit  (he 
characteristic  desired. 

Cnairnller/Llaear  Estimator  lutegreiloa  Results 

The  simple  linear  estimator  was  coupled  with  a  (!cun-Schedu)ing  LQR  Controller 
to  investigate  the  system  perfortTiance.  The  rendezvous  phase  focused  on  remained  to  be 
the  linal  5  km  of  relative  distance  between  tbe  microsaielliie  and  ihe  target,  or  OC-3  in 
Table  2.  Tbe  first  suoulauon  ran  open-loop,  without  raicrosiueUite  thrust  control,  for  200 
minutes.  The  performance  trucked  against  (he  torget  injth  is  provided  m  Figures  32-33 
for  reference.  The  relutive  distance  between  ihe  tnicrosaiellite  and  the  target  truth  is 
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periodic  and  remains  within  20  m  of  the  iniiial  S  ktn  Mparaiion  as  neem  in  the  Ibunh 
8UbpioiofFi|;iire  32. 


Trm  triiiiM^ 


Figure  32.  Open-Loop  LQR^LQE  -  Target  Truih,  Perfect  Initial  Estimate 
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Figure  33  Open-Loop  LQIL/LQE  in  plane  -  Targei  Truih,  Perfeci  Initial 

Esiimaie 

TPe  performance  in  Figure  33  above  doee  not  trad!  to  die  origin  a^  (he  control  input  is 
turned  0^.  The  resulting  open-loop,  or  *'uncontrolle>l.”  perlbrtTiance  is  as  expected. 

The  same  performance  results,  measured  against  die  target  estimate,  are  provided 
in  Figures  34-3 S.  The  relative  distance  between  the  microsatellite  and  the  target  estimate. 
Figure  34,  remains  somewhat  periodic  but  is  more  complex  than  die  perfect  initial 
estimate  case  shown  in  Figure  32.  The  position  of  the  microsatellite  relative  to  the  target 
estimate  captured  in  (he  Sr,  r^Stf  plane  and  shown  in  Figure  35  is  far  ftom  ideal.  I( 
depicts  how  (he  microsatellite  has  a  very  difficult  time  simply  following  the  target 
estimate  in  this  case. 
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Figure  S4.  0(ten-L4Mp  LQR/LQE  •  Target  Esiimate.  Perfect  Iniiial  Eetimate 
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Figure  35.  Open-Loop  LQR/LQE  in  6i\r,S0  pbne  -  Target  Estimate,  Perfect  Initial 

Estimate 

This  result  caused  the  researcher  to  take  a  closer  look  at  the  target  position  error, 
or  the  difference  between  the  truth  and  estimate  The  target  position  error  is  displayed  m 
Figure  3t>  for  this  open-loop  case  The  fairly  periodic  waveform  of  the  error  has 
ampiihide  of  approximately  37  m,  and  period  of  92  minutes.  The  amplitude  was  found  to 
slightly  increase  as  the  simulation  time  was  extended  The  target  state  errors  shown  in 
Figure  36  add  to  the  normally  periodic  open-loop  performance  shown  in  Figure  32 
producing  the  open-loop  estimate  performance  in  Figure  34. 
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Figure  Jf>.  Target  Position  Error,  Perfect  Initial  Estimote 


To  illustrate  ttow  il>e  e«iimaie  error  can  cause  signtficani  difficulties  to  the 
rendezvous  problem,  the  loop  was  closeil  invoking  microsatellkte  control  based  on  the 
current  target  estimate.  The  simulation  again  began  with  a  perfect  initial  estimate  for 
easy  comparison  to  the  open-loop  performance  above  The  results  are  given  in  Figures 
3  7-3  B  below.  The  controller/estimator  developeil  achieves  rendezvous  to  the  target  tn 
204  minutes,  using  9.36  nVs  AV .  It  is  important  to  recall  this  simulation  began  with  the 
microsatellite  trailing  the  target  ftom  only  S  km  and  a  perfect  initial  estimate  of  the  target 
state  Therefore,  results  should  not  be  directly  compared  to  those  by  Tschirhart  outlined 
in  Chapter  II,  but  represent  only  the  final  rendezvous  phase  studied. 
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Figure  37  Closed-Loop  LQItLQE  •  Target  Eeiimaie.  Perfect  Iniciat  Esiimaie 
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Figure  38.  Closed-Loop  LQR^LQE  in  plane  •  Target  Esiimaie,  Perfect  Initial 

Esiimaie 


The  next  la^tk  was  to  include  some  initial  uncertainry  in  the  target  estimate.  TEe 
goal  was  to  demonstrate  if  the  microsatelliie  could  still  rendezvous  to  the  target  starting 
with  a  1  km  target  position  error.  After  running  the  simulation  for  1000  minutes, 
however,  the  minimum  relative  distance  was  never  less  than  39  m.  An  initial  error  of  100 
m  led  to  a  minimum  relative  distance  of  3.6  m,  sdll  greater  than  the  required  1  m 
rendezvous  speciHcation.  Rendezvous  was  achieved,  horn  5  km  out,  in  481  minutes, 
using  12  39  nVs  AV ,  beginning  with  a  10  m  target  position  error.  Although  this  does  not 
represent  an  achievable  initial  error,  it  does  help  characterize  the  integrated 
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comroller/esiimator  perfonnance.  The  results  for  this  final  case  are  shown  in  Figures  39- 
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Figure  39  Closed-Loop  LQR/LQE  •  Target  Eatimaie,  10  m  Initial  Error 
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Figure  40.  Closed-Loop  LQR/LQE  in  Sr,r^66  pUne  -  Target  Esttmaie,  10  m  Initial  Error 

TPe  conclusion  of  integrating  a  linear  filter  with  a  Gain-Scheduled  LQR 
controller  is  (hat  it  is  not  good  enough  to  solve  the  rendezvous  problem  It  only  works  for 
a  nearly  perfect  model  and  initial  guess.  The  main  obstacle  encountered  is  that  a  linear 
filter  cannot  provide  a  good  enough  estimate  given  the  anticipated  iniual  target  error. 

Nonlinear  Filter 

The  on-orbit  Non-(inear  Least  Squares  (NLSl  estimator  developed  by  Foster, 
which  follows  the  seven  step  proces.<  outlined  in  Chapter  IK,  was  used  as  a  starung  point 
for  this  element  of  the  evaluation  The  NLS  filter  was  first  modified  and  then 
characterized  to  understand  the  basic  performance  before  integrating  with  a  controller.  A 
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key  measure  of  fl)!er  pcrfomunce  is  ihe  residue}  term,  or  difTemce  berwecn  meaMired 
and  predicted  observation  values  as  given  in  Ei^uaUon  26 

Ad  open-loop  sunulniion  containing  200  obsen'ations.  separaied  by  one  nuouie 
eukh.  converged  on  an  esumaie  siler  seven  filter  iterauocs.  Tlie  iniiioJ  estinuie  included 
a  ]  km  position  errur.  Tbe  dale  plot  in  Figure  41  displays  the  residual  values  by  iiemiion. 
Subsequent  iierauons  provide  residuals  signiHcandy  better  iban  the  previous,  but  it  can 
be  difljculi  to  ted  the  di^erence  graphically  Tbe  first  iteration  was  removed  in 
generating  Figure  4 1  to  belter  show  the  remaioiog  six.  (t  is  clear  that  iterations  three 
through  seven  ore  significantly  better  than  iterauon  two  Likewise,  removing  iteration 
two  reveals  that  four  through  seven  ore  significantly  better  than  three. 


Figure  41.  NL5  Filter  Residuals  •  Last  Six  Iterations 


73 


The  resldiul  values  fur  Just  ihe  final  iitfraiiun.  in  Fiyure  42,  are  nicely  di  unbilled  urounJ 
aero  aad  cuptured  wtduc  6  m. 


Figure  42.  NLS  Filler  Rasiduats  •  Fuul  Itersuon 


The  resulting  range  residuals  provide  some  confidence  ibe  filter  is  operaiuig  as 
desired.  The  target  posiuun  error  after  the  200  minute  simulaiiun  is  also  of  greui  mteresL 
aa  the  niter-provided  esiimute  is  ibv  mpui  to  the  controller.  The  positloa  error  generated, 
by  iteration,  is  shown  in  Figure  43.  The  initiul  I  km  error  dau  point,  iteration  0,  bas  been 
removed  frara  ibe  chan  for  clarity. 
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Figure  43.  NLS  Target  Posiiioo  Errors 


After  the  Erst  iieraiioo,  the  error  is  reduced  frotn  1  km  to  46  46  m.  When  convergence  is 
achieved,  after  seven  iieraiions,  the  error  has  becrt  decreased  to  1,06  m 

Nununa)  opert-loop  fther  performance  has  been  assessed  to  be  very  good.  When 
the  ft  iter  has  coctverged  on  an  estimate,  both  the  range  resiJoaJs  and  target  position  error 
are  small,  as  desire±  The  nexi  activity  was  to  uiirodiice  microsaieliiic  thrust  cooirol  to 
rende2V0Ui  with  the  target 

Cooiroller/NooUnetr  Estimator  loieyratlon  Elosults 
The  modined  NLS  estimator  characterized  above  was  cuupled  with  a  Gain- 
Scheduling  LQR  controller  to  evaluate  the  potential  for  closed-loop  control  from  a  poor 
mitud  estimate.  The  top-levei  subsystem  arcluteciure  was  developed  based  on  the 
aJgonthn  functions  and  is  shown  in  Figure  44  below. 


75 


Fixture  44.  UN&C  A)gonihfn  AA^biiecniK 


Four  separate  iioie  ueps  were  natded  for  ihe  simulaiion  A  confmDer  ume  stop, 

f  ,  on  esiimaior  ume  siep.  i _ ,  an  obscrvauon  nine  «tp.  ,  and  a 

simuleiiioD  inne,  Table  4  below  descnbas  the  uie  of  each  lime  step.  Tlte  coulroUer 
and  Kiimaior  time  sieps  were  ci^ual  to  eliminate  the  need  for  added  oomplexiiy  to  the 
coniruller  al^^orithm. 


Table  4.  Simulation  Time  Slen  Defuutioiu 


TImeSlv 

&eMier> 

Corrnsllw 

EstinalBtor 

Obaervatton 

SeniJelion 

1  im«  Denwar  inrusT  control  ivdsMa 

Tima  betwear  eslimala  updalas 

Tima  batwaar  maaaurameni  observations 

T  imo  of  otbd  ofopaQallon 
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The  sifnolaiion  ume  u«d  For  ihe  »tafiJ-Qlone  filter  chamcterLmuon  above  svai;  200 


nuDUieiL  The  controller  in  the  iniegraied  sysieni  caocKit  tvaitlOO  nsinuics  for  a  new  target 
e&timate,  however  The  researcher  oeedcd  to  determine  bow  quickly  tbe  filter  could 
provide  a  reasonable  estimate  to  pd.«s  tu  the  cuniroMer  Times  between  the  baseline 
conirulier  (for  perreet  knowledge)  with  a  cuniroller  tune  step  of  60  s  and  the  stand-alone 
filter  lime  of  12000  s  were  considered.  Beginning  with  the  final  5  km  rendezvous  phase. 
OC-3 10  Table  2,  the  integrated  coniroUer’esiinnior  simulation  was  run  for  a  single 
control  loop.  Tbe  resulting  target  position  error  is  plotted  for  various  controller  luoe 
steps  in  Figure  45.  Eliminating  the  first  few  data  points  result  in  a  closer  view  of  the 
remaining  runs,  as  .shown  in  Figure  46. 


Position  &ror  •  Al  CVrta 


0  1030  3COO  3030  dOOO  SCOO  0030  7COO 
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Figure  45.  Simulation  Time  Trade  -  Target  Position  Error 
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Figure  46  Simulatioo  Ttene  Trade  •  Tor^  Position  Error,  Select  Dam 

As  the  coniroiier  controls  based  on  the  rehtiive  posiiioa  between  the  liiT^tei 
estimate  and  tnicrosaielliie.  a  large  position  error  is  dearly  undesirable.  Bused  on  the 
above,  aUowing  die  filter  at  least  900  s  results  in  a  fibrly  good  estimate  Less  time  yields 
a  useless  target  ceircnate. 

Srudyuig  die  relative  distance,  in  addition  to  the  position  error,  exposes  the  fact 
that  the  controller/csiunaior  ume  step  must  be  selected  carefully.  Specifically,  a  dree 
near  1000  s  was  Toimd  lu  be  most  favorable,  for  a  single  control  loop.  This  sirnulation 
mn  reduced  the  relative  distance  from  the  inidal  5.00  km  to  1.04  km  Tune  steps  outside 
the  range  of  700-1400  »  resulted  m  the  microsetcUJte  drifuog  further  away  fiure  the  target 
estimate. 

Having  identified  an  initial  controllcr^esiireaior  lime  step  of  1000  s,  and  an 
observation  time  step  of  10  s,  die  mtcgraied  system  was  simulated  For  multiple  control 
loops  to  an  cITon  to  achieve  final  rendezvous.  Tlte  1000  s  time  step  rumed  out  to  be  too 
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large  for  the  controOer.  however.  The  performance  for  the  first  20  control  loops,  20,000  s 
of  controlled  flight,  is  shown  below  in  Figure  47.  The  sharp  changes  are  a  result  of  going 
too  long  between  thrust  vector  updates  in  the  LQR  control.  A  longer  simulation, 
including  more  control  (oops,  only  results  in  the  rebdve  distance  increasing. 


Figure  47  LQR/N  LS  Performance  for  *  1 000.^ ,  « I  Os 

The  target  position  knowledge  begins  to  significantly  degrade  after  these  first  20 
control  loop  iterations,  a.s  shown  in  Figure  46  below.  This  is  the  cause  for  the  degenerate 
performance  experienced  by  running  the  simubtion  for  longer  periods  of  lime. 
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Figure  4S.  Posiilon  Errof  for  t  s  )OOCl' ,  slO.t 

Ddckmg  off  on  the  Control  Weigtitmg  Mairix,  R ,  was  attempted  to  help  smooth 
oui  die  comrol  vecu^  changes.  Specitkally,  ^  Equations  37  and  38,  was 

increased  lo  lel3,  allowing  less  comrol  thrust  lo  be  applied.  Although  this  did  cause  a 
smoother  rendezvous  approach,  die  filter  was  unable  to  converge  on  an  estimate  during 
the  seventh  control  loop.  Figure  49  shows  the  performance  just  prior  to  filter  failure 
The  performance  for  this  case,  although  still  somewhat  severe  due  to  the  larger  controller 
time  step  of  1000  s,  more  closely  matches  that  of  the  non-gain-scheduled  controller 
assuming  perfect  knowledge  shown  in  Figure  6.  The  large  time  span  between  control 
thrust  updates  has  clear  negative  affect  on  die  abihiy  to  meet  final  rendezvous  critena. 
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Figure  49  LQR/NLS  Improved  PerfomuitKe  for  r  a  )000.« .  *  lO.t 

TlK  other  variable  control  knob,  controllen  estimaior  dme  step,  was  examined 
next.  The  lower  bound  ideniitied  by  the  stmulaiion  lime  trade  study  above,  700  s,  was 
used  to  rerun  the  stmulaiion  This  resulted  in  filter  convergence  failure  ader  only  a  few 
control  neraiions,  keeping  the  same  10  s  observation  time  span  as  above.  When  the 
interval  was  reduced  to  1  &  yielding  700  data  observations,  performance  unproved  The 
smoother  curve  of  Figure  SO,  showing  the  Hrst  12  control  loops,  is  closer  to  that  desired, 
but  still  short  of  rendezvous.  This  140  minute  flight  from  S  km  out  brought  the 
mtcrosaieliiie  to  within  t>7  m  of  the  target. 
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Figure  50  LQR/NLS  I  mpfovc*j  PerfornutKe  for  r  s  700  c .  ‘ 


Although  oui&ide  the  idenii5ed  time  step  range,  a  run  was  conducted  using  300  s 
steps  for  comparison  The  results  for  the  first  12  control  loops  are  illusuwd  in  Figure 
5 1  The  performance  more  closely  corresponds  to  a  smooth  curve,  as  desired,  but  breaks 
down  halfway  to  rendezvous  as  the  estimator  again  fails  lo  converge  on  an  estimate.  The 
above  results  bound  numerous  cases  atiempieil,  none  of  which  yielded  satisfactory 
performance. 
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Figure  51.  LQIWLS  Improved  Performance  for  a  300  c ,  *  Is 

Aliemaiive  CONORS  for  completing  ihe  final  rendezvous  phase  were  explored  lo 
find  a  suIlIcienL  even  if  les«  elegoni,  solution.  One  concept  explored  was  lo  obuin  a 
good  esiimaie  by  following  the  largei  open-loop,  and  then  quickly  controlling  without 
updJtmg  the  estimate.  The  system  would  perform  an  iniual  orbit  determination  using  the 
onboard  sensor  and  then  control  to  chat  solution  without  further  updates  In  this  case,  the 
estimate  containing  a  I  m  initial  error  was  simply  propagated  forward  until  the  catch 
criteria  was  met.  This  scenario  resulted  in  rendezvous  to  the  target  estimate  in  71 
minutes  from  the  S  km  start  The  performance  is  provided  in  Figure  S2.  As  this  case 
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employs  the  some  snull  coniroller  lime  step,  60  s,  as  the  ongino)  LQR.  controller  (hat 
assumed  irtiih,  (he  perfonnance  Is  very  similar  to  that  shown  m  Figure  6. 
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Figure  52.  Blind  Rendezvous  Performance 

Although  the  I  m  relative  distance  criterion,  to  (he  target  eeumate,  was  met,  the 
relative  distance  to  (he  target  truth  is  the  important  measure.  Given  a  small  target 
position  error,  (he  (wo  distances  would  be  very  similar.  In  (his  case,  however,  (he  error 
built  up  over  the  71  minute  simulation  resulted  in  a  12  m  relative  distance  to  (he  utith 
when  the  controller  believed  rendezvous  had  been  achieved 

A  further  complication  with  this  alternative  CX)NOPS  is  (he  lack  of  robustness. 
Once  perturbations  are  introduced,  such  as  air  drag,  it  ^ils  to  achieve  rendezvous  even  to 
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the  tarpei  esiinute.  The  simubdon  was  run,  includjng  J2,  for  200  minutes  resulting  in  the 
performance  stwwn  in  Figure  S3. 


Figure  S3.  Blind  Rendezvous  PerlbrtTunce  with  i2 


In  this  case,  the  target  position  error  becomes  so  large;  there  is  no  hope  of  closure  The 
difference  between  the  estimate  and  truth  continuously  increases  during  the  simulation. 
Figure  54  shows  how  this  error  builds  up  to  over  100  m  al^er^usi  140  minutes. 


8S 


Trm  trriiMi 


Figure  54.  Posiuon  Error  for  Blind  Rendezvous  with  J2 


SumiDtry 

Froni-end  Systems  Engineering  led  (o  a  top-level  system  feasibility,  but  also 
exposed  GN&C  integration  complexity.  Controller  and  estimator  algorithms  have  been 
separately  matureil  previous  work  and  shown  to  contain  a  ^ir  amount  of 
adaptability.  As  an  example,  Gain-Scheduling  was  found  to  significantly  improve  upon 
earlier  Linear  Quadratic  Regulator  performance  Difftculty  arose,  however,  when  the 
development  of  an  integrated  GN&C  subsystem  was  attempted  It  proved  to  be  outside 
the  context  of  established  algorithms  and  required  additional  effort  to  extend 
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V.  Conclusions  god  Recnmmeodatlons 


Cbtpter  Overview 

This  chapter  summarizes  the  leostbiliry  ajialysis  of  a  microsaielliie  rendezvous 
with  non-cooperaiive  targets  and  identiftes  the  technical  challenges  in  doing  so.  The 
tniegraiion  challenges  encountered  in  (he  GN&C  algonthm  assinulaiion  provide  rich 
oppoiTuni(y  for  future  research  activiues. 

C'oocluslons  of  Research 

The  results  of  this  research  underscore  the  difHculry  m  developing  an  iniegroted 
(iN&C  system  for  a  non-cooperaiive  rendezvous.  The  competing  demands  of  (he 
controller  and  estimator  components  prove  (oo  much  for  simple  algorithms  to  overcome. 
The  suite  of  continuous  nonlinear  control,  linearized  dynamics,  and  non-linear 
measurement  conditions  need  to  be  fully  accounted  for  in  the  integrated  solution  The 
NLS  estimator  developed  for  (his  research,  however,  has  no  means  to  include  (he 
continuously  changing  measurements.  This  causes  error  (o  build  up  in  (he  target  estimate 
between  updates,  greatly  hindering  the  process.  Although  a  top-level  system  architecture 
was  developed,  the  technical  complexity  involved  requires  more  sophisticated  methods  to 
solve. 

RcceiQiaendatleiis  for  Future  Research 

As  political  economist  Thomas  Schellmg  has  pointed  out,  "There  is  a  tendency  In 
our  planning  to  confuse  the  unfamiliar  with  (he  improbable.  The  contingency  we  have 
not  considered  looks  strange;  what  looks  strange  Is  ihouglii  improbable;  what  is 
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improbable  ne«d  noi  be  con^dered  sefiousiy.’'  Surprise  is  most  often  not  a  tack  of 
warning,  but  the  result  of  a  tendency  to  di^ia  as  reckle^  wtui  we  consider  improbable 
(Commission  to  Assess  US  National  Secunty  Space  Management  &  Organuaiion,  2001). 
Although  the  conclusion  of  this  research  points  to  a  low  level  of  feasibility,  given  the 
scope,  it  should  not  be  considered  un^miliar  or  improbable. 

As  future  work,  the  control  and  estimation  algorithms  used  need  to  be  further 
tailored  to  incorporate  all  physical  conditions  involved  in  the  rendezvous  problem. 

Taking  everything  into  account  may  reveal  a  feasible  solution.  More  development  dme  is 
needed  to  mature  the  software. 

Although  past  research  indicate  an  impulsive  thrust  solution  does  not  work,  even 
given  perfect  knowledge  of  the  target,  it  seems  additional  study  may  reveal  a  viable 
produci.  As  it  was  very  difficult  to  integrate  a  continuous  thrust  controller  with  an 
estimator  in  this  work,  the  impulsive  controller  should  be  revisited.  Specidcally, 
coupling  a  CW  controller  with  a  modified  NLS  Riter  could  present  interesiing  results. 
Optimal  sequenoa)  processing  with  noise  statistics  assumed  should  be  investigated 

Once  a  technical  solution  is  found,  the  systems  engineering  work  should  be 
expanded  to  include  a  more  formal  evaluation.  A  full  system  architecture  should  first  be 
developed  followed  by  more  rigorous  evaluation  techniques.  Executable  models  could  be 
developed  upon  which  to  base  a  quantitative  assessment  of  system  feasibility.  A  detailed 
ON&C  subsystem  feasibility  study  should  be  conducted  to  explore  hardware/software 
integration  challenges,  followed  by  a  full  system  study  to  look  at  subsystem  integration 
issues. 
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AppeiidU  A  -  Main  LQR  Code 


% 

%  MAfN  LQR  CODE 
% 

%  Ailen  Toso  -  7  Feb  04 


% 


%  Orbiia)  Rendezvous  Wiih  a  Non-Cooperanve  Target 
% 

%  Onginal  code  provided  by  Troy  Tschirrtan  (4  Mar  03) 

%  Code  tnodiEed  (o  incorporate  Esdmaior  and  focus  on  final  5  km  approach 


% 

% 


% 

%  This  program  uses  (he  following  function  files  which  must  be  on  (he  current  path' 
% 


%  atmosphere.m 
%  Calcintt.m 
%  D0_Pl0(8 
%  LQR_Rend 
%  propagate.m 
%  posvel.m 
%  kjk2pqw.m 
%  pqw2ijk.m 
%  n22pqw.m 
%  rv2coe.m 
%  coe2rv'.m 
% 

V  ■  ■  ■  >  ■  i 


calculate  aimosphenc  density  ai  the  given  atiiiude 
calculate  the  Initial  conditions  for  the  run 
plot  (he  results 

accomplish  Iqr  rendezvous  manuevee 
propagator 

sei  up  the  ditTerenilal  equation  for  the  propagator 
transform  r,v  from  ijk  frame  to  pqw  frame 
transform  r,v  from  pqw  frame  to  ijk  frame 
translbrm  r,v  from  rtz  frame  to  pqw  frame 
calculate  coe  for  the  given  r,v 
calculate  r,v  for  the  given  coe 


% . I  ■  .  — - 

%  Clear  Variables  and  Set  Format  Options 

% — . . 


clear 

format  long  g 
format  compact 

%  Print  a  banner  to  separate  results 

%  Starr  the  timer  (used  at  the  end  to  deiermine  how  long  the  run  took) 
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tic 


%  Sei  Selectable  Variable  Values 


% . . 

%  Trve  target's  actual  initial  COEs 
•4 . . 


coe_t8t_act(  1)  -  6  772RBB9l22MB40e«*003; 
coe_tet_act(2)  -9  BB77l3549B259l3e-004; 
coe_t8t_9Ct(  3 )- 0  797363  B64B5  B27; 
coe_t8t_act(4)  -  0  9O75799OO7B3B0; 
coe_t8t_9Ct(5)  -  I  5IB437609B069I; 
coe_tgt_act(6)  -  5  59054044657763; 
coe_tgi_eci(7)  -0  0; 


%kin  a 

%  dimensionless 


%  radians 
%  radians 
%  radians 
%  radians 
%  seconds 


nu 

I 

cap_omega 
small_omega 
time  since  perigee 


•4 . . 

%  Set  initial  micro  offset  from  the  target 
. . 

deir  -  0;  %  kilometers  (delta_r) 

dist  -  -5 ;  %  kilometers  arclengih  (ro”delta_ineta) 

delz  •  0;  %  kilometers  (delia_z) 

% . . 

%  Set  the  acceptable  relative  distance  and  velocity  for  a  successful  rendezvous 
% . . 


catcMis  -  0  00 1 ;  %  kilomeieis  ( 1  m)  -  Phase  2  (0/L),  Phase  3  (OL) 

catchve)  *  0.00001;  %  kilometers/second  ( I  cm/ si  •  Phase  2 


timestep  -  1000;  %  used  onty  to  determine  delta_thnjst  in  LQR,  based  on  filter 

%  (seconds) 


•4 . . 

%  Specify  values  for  ihe  state  weighting  matrix,  Q 
%  and  the  control  weighting  matrix,  R 

%  Note:  Q_mag  increases  «>  faster  movement  from  initial  to  desired  states 
%  R_mag  increases  ->  lower  control  usage 
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% 


Q_mag  -  I ; 

R  fin  -  Iel3; 


% . . 

%  Sei  (temirbaiion  (J2)  opiuns  (Note!  "1"  -  option  selected;  *'0*'  -  option  noi  selected) 
% . . - . 


pen-0; 


% . . 

%  Sei  drag  optN>ns  ond  valuer  (Note;  "1"  -  option  selected;  'O*'  ~  option  not  selected) 
•4 . . 


dragtgt-O; 

cdjgt  -  2.2;  %  Drag  coefTicleni  of  the  target 

ajgt  -3.5*1 .2/(  1 000*2);  %  Area  of  die  target  ( l:m*2 ) 

mjgt  -  725;  %  Mass  of  the  target 

cdamtgt  -  dragtgt "  (cd_tgt "  a_tgt)  /  m_tgt;  %  Calculate  the  target's  cdam  value 

•4 . . 


dragmic  -  0; 

cd_mic  -  3;  %  Drag  coefficient  of  the  micro 

a_mic  -15/(1 000*2 );  %  Area  of  the  micro  (km*2) 

m_mic  -  1 00;  %  Mass  of  the  micro 

cdanmic  -  dragmic  •  (cd_mic "  a_m»c)  /  m_mic;  %  Cakubte  the  micro's  cdam  value 


. . 

%  Set  plot  options  (Note;  'r  -  option  selected;  ‘O "  -  option  not  selected) 
% . . 


prdijk  -  0;  %  Plot  relative  distance  in  the  inertial  (ijk)  fi’ame 

prdrtz  -  1;  %  Plot  relative  distance  in  the  relative  (nz)  fiame 

prdroto  -  I ;  %  Plot  relative  distance  in  the  relative  plane  (delta_r,  ro’'de!ta_theta) 


% - 

%  Initialize  variable  values 

% . - . . 


delta  V  accum  -  0;  %  initialize  delta-V  to  zero 


%  Calculate  Initial  Values 
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%  I  Targei's  initial  position  and  velocity 
%  2  Micro's  initial  position  and  velocity 

Calcinit 


%  Accomplish  Linear  Quadratic  Regulator  Rendezvous 


LQR_Rend 


. 

%  Print  Output  Values 


delta_r  -  dels 
ro_delta_theta  ~  dist 
delta_2  -  delz 

pert 

dragtgi 

dragmic 

P' 

time  step 

Q_mae 

R_ma8 

delta_v_accum 
final_vel  -  vel_iiow 
Rnal  dist '  dist  now 


%  Draw  Desired  Plots 


Do  Plots 


%  End  of  Program 


run  time  -  toe 
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App^ndlY  B  -  Lio«or  Estimator  Code 


%■■  ■■ 

% 

%  STAND  ALONE  LQE_LSIM  ESTIMATOR  (Linear  Estimator) 
% 

%  Ailen  Toso  -  2  Dec  03 
% 

%  Orbiial  Rendezvous  With  a  Non-Cooperaiive  Target 
%  Uses  Range.  Az.  El  to  obtain  an  estimate  for  the  target  state 
% 

%  (L,P.E)-LQEtA.G,C,Q,R) 

% 

%  A  -  Plant  (System) 

%  G  ~  something  small 

%  C  ~  Observation  Geometry  Matrix  -  Range  Vector  Range,  Az.  El) 
%  Q  ~  Strength  of  Process  Noise,  set  to  zero 
%  R  '  Strength  of  Measurement  Noise 
% 

%  [y.x]  -  LSJM(A.B.C.D,u,l.xjet') 

% 

%  [y_hat,  x_hatj  -  LSIM(A_ob,L,C,D,y,i,x_tgt_hat’) 

% 

%  Initial  input  % 
clear 

formal  long  g 
formal  compact 

xjgt  -  [0  1  0  0  0  OJ 

r_tgi_reKl)-x_tgU2>; 
r_tgi_rel(2)-x_tgf(4); 

v_tgt_reKl)-x_tgi(I>; 
v_tgt_rel(2)-x_tgt(3>; 


x_tgt_hat  •  x_tgt; 
p08_error  -  .70710678; 


%  detia  r  only 


%  amounts  to  I  km  position  error 
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f_igt_hai_re((I )  -  x_igi_fwt(2)  -  pos_error; 
r_igt_hai_re((2)  -  x_igi_fwt(4)  -  pos_error; 
r_igt_hai_re((3)  -  x_tgi_hat(6) 

v_igi_hai_re((l )  -  x_tgi_fwi(l ); 
v_igt_hai_re((2)  -  x_tgi_lwi(3); 
v_igt_hai_re((3)  -  x_tgi_hat(5); 

x_iet_hai(2)  -  r_igi_hai_rel(  I ); 
x_iet_hai(4)  -  r_igi_hai_rel(2); 
x_iet_hai(6)  ^  r_iel_hai_rel(3) 


%  Calc  A 


U  -  39SA0 1 :  %  km'3/8ec'2 

coejgl_aci(  1)  -  6  772  BB  89 1 22MB40e-MM3;  %  a,  km 

coe_ref-  coe_tgt_dct;  %  Stan  with  target's  Initial  COEs 

n-sqrt(U/coe_ref(iy‘3); 

A-(0  3"n'2  2"n  0  0  0; 

10  0  0  0  0; 

-2*n  0  0  0  0  0; 

0  0  10  0  0; 

0  0  0  0  0  -n'2; 

0  0  0  0  I  Oj; 


%  Calc  C 


%  Calc  the  norm  of  target  uaie,  x_tgi,  in  nz,  in  km 

%  Relative  position  vector  (3x1)  (range)  in  nz  coordinates 
%  to  die  target  satellite. 
fange_veciDr  •  r_tgt_hat_fei; 

%  Magnitude  of  the  range  vector  in  nz,  in  km 
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range  -  norm(raiige_vecior); 


%  C  ~  Ob^rvacton  Geometry  Matrix  -  Range  Vector  Range,  Az.  EO 
%  Intiialize  C  (o  zeros  Hrsi  (hen  build  up  needed  components 
C  -  zeros(3.6); 

%C(l.l).C(2.l).Q3.l)-0 

%C(U).C(2^).C(3^)i 

C(  I  ^ )  -  fange_vecior(  1 )  /  range; 

Q2^)-  (-fange_veciorl2yfange_veciorll)’'2)/(l  • 
(fange_vectDr(2yrange_vectDr(  1  ))''2); 

C3_2_iop  ^  (-range_vecior(Ii*range_vecior(3))  /  ((fange_vecior(  1)^2  • 
Kinge_vecior(2)''2)’<3/2)); 

C3_2_bonom  -  I  •  (range_vector(3>'2i  /  (range_vector(iy'2  •  range_vector(2r2); 
C(3^)  -  C3_2_top  /  C3_2_bonom; 

%C(U).C(2^).C(3^)-0 

%QI.4).Q2.4).Q3.4)i 
C(l.4)  -  fange_vecior(2)  /  range; 

Q2.4)- (I/range_vecior(Iii /(I  *  (range_veciorl2Vrange_vecior(l)r2); 

C3_4_iop  -  (-range_vector(2i"range_veciorl3))  /  ((range_vecior(  1)^2  • 
Kinge_vecior(2)''2)’<3/2)); 

C3_4_bor!om-  I  •*  (range_vector(3>'2i  /(range_vector(iy'2  •  range_vector(2r2); 
C(3.4)  -  C3_4_top  /  C3_4_bonom; 

%C(l4).C(24).Q34)-0 

%QI.6).Q2.6).Q3.6)i 
Q  1.6)  -  range_vecior(3)  /  range; 

C(2.6)-0; 

C3_6_iop  ^  I  /  ((range_vector(l>'2  ••  r4nge_vector(2>'2r{ly2)); 

C3_6_bonD(n-  I  ••  (range_vector(3>'2)  /(range_vecior(I>^2  •  range_vector(2r2); 
C(3.6)  -  C3_6_top  /  C3_6_bonom; 


%  Calc  L,  esiimaior  gain  matrix 


%  [L,P.E)-LQE1A.G,C,Q,R) 
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%  A  -  Plant  (System) 

%  G  *  vmneihing  small 

G-[l  1  1; 

I  I  h 
I  I  I; 

I  I  h 
I  I  h 
I  I  1): 

%  Form  Q,  the  [irocess  noise  nutnx  nonzero  •  two  ordefs  of  nugnicude  smaller  than  R. 
allow  R  (o  dominate 

%Q(1,1)  -0.0002*2;  -  Reference  Foster's  obser  m 

Q-2eros(3,3>; 

Q(l.l) -0.00000004; 

Q<2,2)  -  0.00000004; 

0(3,3)  -  0.00000004; 

%  R  •  Strength  of  Measurement  Noise 
%  Form  R.  the  jnstrumenial  covariance  matrix 

%  R(  l.l )  -  0.002*2;  Instftjmeniaiion  sigma  squared  (  2  meters  -  0.002  km) 

R-2eros(2^); 

R(  l.l  >-0.000004; 

R(2.2)- 0.000004; 

R(3.3)- 0.000004; 


[L,P,E)  -  lqe(A,G.C.Q,R); 


%Can  LSIM 

%— .  i  ■  i  i  . 

B-JO; 

0; 

0; 

0; 

0; 

oj; 
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D-0; 

A_ob  -  A-L-C; 
I  -  (O' 1:10000]'; 


[yAj  -  lsun(A,B.C.D,u,UJ8’'); 

[y_hjt,  x_haij  -  lsim(A_ob,L,C,D,y,t,x_tgi_hac*); 


V, - 

%  Ploi  resulci 

% . 


%  Parse  x 

rjet_ploi(;,l)  -  x(:,2>; 
rjet_ploi{;,2)  -  x(:,4>; 
rjet_ploi(;,3)  -  x(:,6]; 

vjet_p(oi(:,l)  •  x(:,I>; 
vj8i_p(o«;,2)  -  x(:,3>; 
vjglj5loi(;,3)  -  x(:,5>; 

%  Parse  x_hat 

r_igt_hai_plot(:,) )  -  x_hai('^); 
r_igt_hai_plot(:i)  -x_hai('.4); 
r_igt_hai_plot(:,3)  -  x_hai(  .6); 

v_igt_hai_plot( 1 )  -  x_hai( .  I ); 
v_igt_hai_plot(:^)  -x_haU'^); 
v_igt_hai_plot(:,3)  -x_hai(  4); 

fieiJre(l);  clf; 

pkpi(i/_igi_hai_plou.f_t8«_pfoi.'**'>;8'’'^  o";  ylabel(  Target  Position  Error  (km)'); 
xlabelCTime  (minutes)'); 

fieiire(2);  clf; 

p(oi(i,v_igi_hai_plou.v_l8ij>krt, '--');  grid  on;  ylabel(Targei  Velocity  Error  (km)'); 
xlabel(Time  (minutes)'); 
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%  Check  Qrgei  posiiion  and  velociiy  errors 
V,— - — - - - 

x_igi_nexi  -  x(  10000.:) 

r_igi_rel(  1)  -  x_igi_nexi(2); 
r_igi_reK2)  -  x_igi_nexi(4); 
r_igi_reK3)  -  x_igi_nexi(6) 

v_igt_rel(  1 )  -  x_igl_nexi(  1 ); 
v_igt_rel(2)  ^  x_iet_nexi(3); 
v_igt_rel(3)  ^  x_iel_nexi(5) 


x_igt_hai_tKxi  -  x_hai(  10000,0 

r_igi_hai_re((I>  -  x_igi_hai_next(2>; 
r_igt_hai_re((2)  -  x_igi_twi_nexi(4); 
r_igt_hai_re((3)  -  x_igi_hai_nexU6> 

v_igt_hai_re((l )  -  x_tgi_h4i_nexUl); 
v_igt_hai_re((2)  -  x_tgi_hat_nexU3); 
v_igt_hai_re((3)  -  x_tgi_hai_nexi(5> 


igt_posn_ertor  -  r_tgt_rer  -  r_tgi_hai_re(' 
igt_posn_ertor  •  nonTjtgi_posn_error) 


igt_vel_ertor  •  v_igt_rer  -  v_tgi_hai_re(' 
igi_vel_ertor  -  twrm(igi_vel_error) 


% - 

%  End  of  Program 

- - 
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lAMSTIU^ 

The  feasibiliiy  of  using  a  niicrosaiclliie  lo  accomplish  an  orbimJ  rendezvous  with  a 
QOQCooperauve  target  was  evalualed.  This  study  focused  on  idenulying  and  lunber 
explonog  ibe  lechcical  challenges  involved  in  achieving  a  ooocooperaiive  reodezv'ous. 
A  system  engioeering  analysis  and  review  of  past  research  quickly  led  lo  a 
conceairetion  on  die  guidance,  navigalion,  and  control  elements  ofibe  microsaielhie 
operation.  The  iniegrauon  of  control  and  orbii  deierminaticn  algonibms  was 
invesiigaied.  A  simple  yet  robust  suluiion  could  not  be  found  lo  meet  reasonable 
rendezvous  ctiicrio.  using  essentially  oR'-ihe-shelf  lechnology  and  algorithms.  Syaem 
feasibility  has  been  aiisussed  lo  have  a  low  probability  in  the  very  near  term. 
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